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Preface 


Ground-based /terrestrial radar interferometry (GBRI) is a scientific topic 
of increasing interest in recent years. The GBRI is used in several field as 
remote sensing technique for monitoring natural environment (landslides, 
glacier and mines) or infrastructures (bridges, towers). 

These sensors provide the displacement of targets by measuring the phase 
difference between sending and receiving radar signal. If the acquisition rate 
is enough the GBRI can provide the natural frequency, e.g. by calculating 
the Fourier transform of displacement. 

The research activity, presented in this work, concerns design and devel- 
opment of some advanced GBRI systems. These systems are related to the 
following issues: detection of displacement vector, Multiple Input Multiple 
Output (MIMO) and radars with 3D capability. 

The conventional GBRI measures only the component of displacement 
along range direction. A GBRI operating in monostatic and bistatic modal- 
ity is presented in this work. The sensor detects the first component of 
displacement as the conventional GBRI (monostatic) and an additional com- 
ponent through a transponder (bistatic). 

The radar has been successfully tested in controlled environment using 
a basic transponder (two antennas and an amplifier). The transponder has 
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been improved in order to increase the gain of the amplifier and to solve some 
issues of the basic version. Finally, the system is used in real application for 
measuring the natural axis of a telecommunication tower. 

The most advanced GRBI system can measure the direction of arrival 
of scattered signal by exploiting the movement of the antenna on an axis 
(Ground Based Syntetic Aperture Radar - GBSAR). The step between two 
position on the axis has to be smaller than a quarter of wavelength. 

The emerging Multiple Input Multiple Output (MIMO) technique can 
be used to reduce the mechanical movement parts and the problems related 
to these. Also for MIMO radar the spacing between two closer phase center 
has to be smaller than a quarter of wavelength for the Shannon theorem. 

In this work a Compressive Sensing (CS) MIMO radar is described. In- 
deed the CS is a technique able to reconstruct signal without the constrain 
of Shannon theorem. The signal has to be sparse and randomly sampled in 
order to use the CS. 

The CS technique can be applied for increase the scan-length of a MIMO 
system of 40% + 50%. Therefore, by using the same number of antennas, 
the CS allows to increase the angular resolution of a MIMO radar. 

A prototype of interferometric CS MIMO radar has been developed and 
tested on some bridges. The results were compared with a conventional 
GBRI with a good agreement. The CS MIMO radar was able to discriminate 
the left-right movement of bridges. Unfortunately the repetition rate of this 
prototype was not enough to retrieve the spectra of natural frequency. 

Since the movement is along a single axis the obtained radar image does 
not have angular resolution in the plane orthogonal to the scan axis. In 
other words, if the radar head scans along the x-axis the radar image cannot 
have resolution in elevation angle. This is not a serious problem when the 
scenario is a slope, where the elevation (z-axis) can be reasonably considered 
an unambiguous function of the (x,y) position. Unfortunately there are cases 
where the geometry of the structure under test is much more complex, i.e in 
urban environment. 

In this thesis two radar systems with three dimensional resolution are re- 
ported. This two systems synthesize the two technique previously described. 
Indeed the first sensor uses the bistatic principle by exploiting the move- 
ment of an additional antenna in vertical axis for obtaining the resolution in 
elevation. The second system exploits the movement on an horizontal axis 
of the CS MIMO with phase center positioned on a vertical axis. 

In order to test the capability, the two radars were located in a urban 
scenario in front of a 7-storey building. Both systems were able to provide 
a 3D image of the building. 


State of the art of Ground-based radar interferometry 


Ground-based /terrestrial radar interferometry (GBRI) is a scientific topic of 
increasing interest in recent years. The GBRI is a popular remote sensing 
technique for monitoring landslides, mines, bridges, towers, dams, and other 
civil infrastructures. 


The radar interferometry started as space technology. In the nineties, the 
satellites ERS-1, JERS-1, RADARSAT-1, and ERS-2 [1], [2], made use of the 
phase information of radar images for detecting ground changes. Spaceborne 
Synthetic Aperture Radar (SAR) systems, operating from an orbit at 800 km 
altitude, provide impressive interferograms due to the ground displacement 
after an earthquake [3], [4]. These extraordinary spaceborne developments 
had an early follow-up in analogue ground-based radar systems. Pioneering 
works anticipated some ideas. In 1997 [5], Tarchi et al. detected by radar 
the structural change of a beam in a large anechoic chamber. In 1999 [6], the 
same bulk laboratory equipment was used to detect the changes of a dam. In 
2000, Pieraccini et al. [7] measured the deformations of a pedestrian bridge 
during its static test with a portable radar. In 2003 [8], a landslide was 
monitored by an interferometric Ground-based Synthetic Aperture Radar 
(GBSAR) for the first time. 


In February 2000, a specifically modified radar system with one trans- 
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mitter and two receivers (the second at the end of a 60-meter mast that 
extended from the payload) flew on board the Space Shuttle Endeavour dur- 
ing the 11 day STS-99 mission [9]. This equipment was able to generate a 
high-resolution digital topographic database. Inspired by this idea, a GB- 
SAR has been modified for acquiring data with different baselines to obtain 
the Digital Elevation Model (DEM) of its field of view [10]. 

The strength of GBRI is its complementarity with spaceborne SAR [11], 
[12], [13], while the latter is a powerful and effective remote sensing technique 
for surveying large areas (many km?) at long term (return time in the order 
of several weeks or months), GBRI is suitable for monitoring small areas 
(from a single building to a slope) at short term (sampling time up to a few 
of minutes) as sketched in Figure 1.1. 


Minutes Hours Days Months Years 


Figure 1.1: Complementarity of GBRI with spaceborne SAR. 


1.1 GBRI in the Scientific Literature 


Since the pioneering works in the nineties, the scientific literature of GBRI 
has greatly grown involving several academic groups, private companies, 
and government agencies. Figure 1.2 shows the number of publications from 
1997 [5] to December 2018. The average bibliographic production of the last 
four year is 42 publications per year. 

In the following section the publications are classified in terms of opera- 
tive modality (design and operative bands) and in terms of its applications 
(landslide, bridge, building, mine,...). Finally some of the recent develop- 
ments are briefly discussed: detection of displacement vector, Multiple Input 
Multiple Output and three-dimensional image. 
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Figure 1.2: Papers published until December 2018 on GBRI [14]. 


Radars can detect a range by measuring the time-of-flight between the 
transmitted and received electromagnetic wave. The direction of arrival 
(DOA) cannot be detected (briefly “no-DOA”), [15], [16], [17], [18] or can be 
detected by the rotation of the radar head (“Rotary”) without any synthetic 
aperture [19], [20]. The most advanced radar systems are able to detect the 
direction of arrival through synthetic aperture (SAR: Synthetic Aperture 
Radar). The main advantage of SAR compared with equal aperture (physical 
or scanned) is that the angular resolution is twice good [21]. The simplest 
Ground-based Synthetic Aperture is based on the movement of a radar head 
along a linear mechanical guide (“linear SAR”) [7], [8]. Alternatively, it 
can be based on a circular movement (e.g., along an arc) [22], [23], [24]. 
This kind of radar is called “C-SAR”. Finally, the synthetic aperture can be 
made with multiple transmitting antennas and multiple receiving antennas 
(MIMO: Multiple Input Multiple Output) [25], [26]. A less common way 
to realize a synthetic aperture is to use the movement of a helicoidally slot 
[27]. Figure 1.3 shows pictures of some prototypes that exemplify these six 
operative modalities. 

As shown in the histogram in Figure 1.4, the majority of the papers 
refer to linear SAR (58%) [7]. About 25% of the papers discuss develop- 
ments or applications of radar interferometers without angular resolution 
(“real aperture”). About 8% of the publications refer to radar systems that 
achieve angular resolution just by rotating the antenna. These three opera- 
tive modalities cover 91% of the publications. 

The operative band is a key parameter of a radar. Generally speaking, by 
increasing the operative frequency, the angular resolution improves and the 
displacement resolution increases. At the same time, the temporal coherence 
of the targets decreases. Therefore, a balance is necessary. Figure 1.5 shows 
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(a) no-DOA 


(d) MIMO (e) ArcSAR (f) moving slot 


Figure 1.3: Examples of ground-based radar interferometers: (a) no-DOA 
(after [28]); (b) Rotary (after [29]); (c) linear SAR (after [28])); (d) MIMO 
(after [30]); (e) ArcSAR (after [48]); (f) moving slot (after [27]). 


the histogram of the publications divided on the basis of the chosen operative 
frequency of the radar. About 58% of publications are relative to equipment 
in Ku-band (12 - 18 GHz). Some radar systems operate in C-band (4 - 
8 GHz), and X-band (8 - 12 GHz). The Ku-band is a good trade-off between 
size of the equipment and temporal coherence of the targets, even if the C- 
band is probably preferable for monitoring snow covers [31], [32], [33]. The 
components cost can be another issue that can make the difference in the 
selection of the band. Until today, components in Ka and W were more 
expensive, but with the large diffusion of automotive radar in W-band, their 
cost is rapidly decreasing [34]. 


GBRI is an applicative topic. Figure 1.6 shows the whole the publications 
grouped based on the specific application. Obviously, a single publication 
can cover more than one application, and there are publications that do not 
refer to any specific application. 


Landslide monitoring is doubtless the most popular application of GBRI 
(see Figure 1.7). Sinkhole monitoring [35] has been included in this category. 
Over 40% of the publications are developments or case-studies about mon- 
itoring systems of landslides. The first in-field test dates back to 2003 [8]. 
Landslide monitoring is an application of great social and scientific impact 
and currently many slopes are routinely monitored by radar interferome- 
12 
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% of publications 


no DOA Rotary Linear SAR C-SAR MIMO moving slot 


Figure 1.4: Papers published until December 2018 classified on the basis of 
the operative modality [14]. 


% of publications 


20 


l a n 
o | MI — imi | — 


L1-2GHz S 2-4 GHz C 4-8GHz X 8-12GHz Ku12-18GHz K18-27GHz Ka 27-40 GHz V 40-75 GHz W 75-110 GHz 


Figure 1.5: Papers published until December 2018 classified on the basis of 
the operative band [14]. 
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Figure 1.6: Papers published until December 2018 classified on the basis of 
their applications [14]. 


ters [36], [37]. 

Testing or monitoring of bridges is another key application (see Figure 
1.7). The radar is a fast, reliable and accurate equipment for remotely 
sensing the static or dynamic behaviour of a bridge. Pieraccini et al. were 
the first (2000) that demonstrated [7] this possibility. Overall about 30% of 
publications refer to bridge monitoring. 

Building is the third class of applications in order of popularity (see 
Figure 1.7). This class is less homogeneous if we considered radar systems 
that monitor construction works [41], villages built over landslide [42] or 
on terrain subject to subsidence [43]. Probably, this class of application is 
overestimated because many papers deal with the tests of equipment rather 
than their effective applications. 

On the contrary, the class of the open pit mining is probably underes- 
timated. Indeed, hundreds of radars are currently in operation as safety 
equipment in mining areas (see Figure 1.7), but the number of publications 
on monitoring activity for worker safety in mining areas is rather modest 
(about 15%). This is probably because mine safety is a rather confidential 
issue, and the owners hardly agree to publish their radar data. The first 
publication about it dates back to 2000 [19]. 

Monitoring of towers/chimneys/electric poles is another emerging appli- 
cation field of GBRI (see Figure 1.8). Generally speaking, all nearly one- 
dimensional structures, like bridges and towers, are particularly suitable for 
14 
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Figure 1.7: Examples of applications of GBRI: (a) landslide monitoring (after 
[38]), displacement map and picture projected on the digital elevation model 
of the slope; (b) monitoring of a viaduct in Florence, Italy; (c) displacement 
map projected on the picture of the monitored building (after [39]); (d) 
mines monitoring, radar deployed in Thompson Creek Mine, Idaho, USA 
(after [40]). 


being monitored/tested by radar, as it does not require cross-range resolu- 
tion. Therefore, if radar operates in no-DOA modality the structure can 
be sampled at over 50 - 100 Hz [15], [16], and its dynamic characteristics 
(natural frequencies, modal shape, transient response) can be detected. 

Glaciers are fundamental “proxies” for standing the effects global warm- 
ing on a local scale [48]. So, scientists have great interest in using GBRI for 
monitoring the change of glaciers (see Figure 1.8). The first case study dates 
back to 2007 [49], but aftermath many other studies on different glaciers have 
been performed (see for example [13], [50], [51]). In our bibliographic survey 
about 13% of the papers refer to glacier monitoring. 

Another field of application of great scientific and social interest is the 
monitoring of cultural heritage (see Figure 1.8). Radar has been used for 
monitoring the lending Tower of Pisa, Italy [16]; the towers of San Gimignano, 
Italy [52], the Bell-tower of Giotto in Florence [53], the Michelangelo’s David 
at Accademia Museum in Florence [46], the Baptistery of Florence, Italy [54], 
the ruins of Roman Forum of Rome, Italy [55], the retaining walls of the his- 
torical town of Volterra, Tuscany [56], and other sites [57], [58]. Overall the 
publications the refer to Cultural Heritage represent about 10% of the total. 
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(b) glaciers 


(c) monuments (d) volcanoes 


Figure 1.8: Examples of GBRI applications: (a) towers (after [44]); (b) 
glaciers (after [45]); (c) monuments (after [46]); (d) volcanoes, permanent 
installation at the Stromboli volcano, Italy (after [47]). 


Volcanoes are both a fascinating phenomenon and an impending risk for 
people living around their flanks. Therefore, although the number of active 
volcanoes is rather limited, over 8% of the scientific papers refer to volcanoes. 
In particular, Stromboli volcano, Italy, (see Figure 1.8) is one of the most 
studied volcanoes in the world [59], [60], [61], [62], [63] and it has been 
continuously monitored by terrestrial radar since 2008 [64]. Also, Soufrière 
Hills in the Montserrat island [65] and Hakone volcano in Japan [66] have 
been monitored by terrestrial radar. 

Digital Elevation Model (DEM) production is another interesting GBRI 
application (see Figure 1.9). Twelve papers (6% of total) have been published 
on this subject. The main finding is that this technique is effective even 
by a terrestrial radar, but its accuracy is sensibly worse than traditional 
equipment (photogrammetry or Terrestrial Laser Scanner) [10]. 

Dam-monitoring has been the first in-field application of GBRI: it dates 
back to 1999 [6] (see Figure 1.9). The concrete surface of a dam is almost 
an ideal target for radar interferometry; it is coherent, homogeneous enough 
but not specular. Nevertheless, the papers related to this application are 
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(c) snow cover (d) trees 


Figure 1.9: Examples of GBRI applications: (a) terrain mapping, digital 
elevation model of a slope obtained through GBSAR (after [67]) (b) dams 
monitoring (courtesy of IDS georadar [28]); (c) snow cover monitoring, map 
of snow water equivalent change obtained by GBSAR (after [68]); (d) trees 
monitoring by GBSAR (after [69]). 


only 5% of total, but dams have already sophisticated monitoring systems 
with different sensors [70], so there is no demand for novel sensors; and dam 
safety, as well as mines, is a rather confidential issue, and the owners hardly 
agree to publish their radar data. 

Martinez-Vazquez in 2005 [31] proposed the GBRI for assessing snow 
cover change (see Figure 1.9), as it has already been done by satellite [71]. 
The technique’s results are effective with dry snow, while wet snow can be 
a problem [32]. Very interesting results have been obtained using the radar 
for mapping the avalanches on steep slopes [33]. 

Monitoring of vegetation or trees by GBSAR (see Figure 1.9) [69] is again 
a marginal application (2% of the total). 

Finally, it is worth noting two very peculiar applications. On January 
2012, the Italian cruise ship Costa Concordia ran aground and overturned 
after striking an underwater rock close to the Isola del Giglio, Tuscany. The 
wreck came to rest on a rock ledge several months before its translation and 
decommissioning. In this period, an interferometric radar was installed over 
a rock of the island for monitoring possible displacements of the wreck [72]. 
The radar detection of the displacement of a pier pulled by a tugboat is the 
second noteworthy application [73]. 

These are the applications that we can find in the scientific literature, but 
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they surely do not exhaust all the possibilities of GBRI. Just as an exam- 
ple of a possible interesting development, we can refer to some preliminary 
works on the use of radar for monitoring wind turbines [74], [75]. The ro- 
tation of the blades gives a peculiar doppler signature (named “sinogram” ); 
Li et al. [76] had the idea to apply the ISAR (Inverse Synthetic Aperture 
Radar) techniques for focusing the sinograms. In this way they obtained 
high-resolution radar images of the blades. A further step could be process- 
ing interferograms between different rotations for detecting the deformation 
of the blades under stress. The same technique could be used for monitoring 
any other turbines like steam turbines and turbofans of airplanes. 


1.2 Multiple Input Multiple Output (MIMO) 


The linear GBSAR as well as any rotary radar and C-SAR need the physi- 
cal movement of a radar head. The mechanical moving parts are the most 
prone to failures, especially in systems that operate outdoors for long peri- 
ods. This is the case with the monitoring of mines and slopes. The phased 
array radar [77] is a possible solution. This technology is well-known and 
widely used in the military field. Its problem is the cost, which is too high 
for civil applications. An emerging paradigm in the radar field is the MIMO 
(Multiple Input Multiple Output). The idea is to use a set of Nrx trans- 
mitting antennas and a separate set of Nrx receiving antennas suitably 
positioned in space (see Figure 1.10). 


Virtual antennas 


Tx Tx, RxRXRXRx, 


n 2 


Figure 1.10: MIMO working principle, Nrx transmitting antennas with 
Nrx receiving antennas are equivalent to Nrx x Nrx “virtual” antennas. 


For each transmitting antenna, all the receiving antennas are working (si- 
multaneously or sequentially), so the number of independent measurements 
is given by N = Nrx x Nrx. The resulting system is a radar equivalent to 
a phased array that needs much less transmitting/receiving channels. 

The first to propose this approach for a GBSAR was Hong in 2010 [78]. 
Tarchi et al. published in 2013 [25] a notable implementation of MIMO- 
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GBSAR named MELISSA; 16 transmitting antennas operating with 16 re- 
ceiving antennas, an equivalent aperture of 2.56 m and an acquisition time 
of 3.6 ms. This very high speed of acquisition opens new interesting perspec- 
tives to this kind of radar, especially in the dynamic monitoring of buildings, 
towers and bridges. The MIMO approach, a radar like MELISSA, is again an 
equipment rather complex and expensive, that hardly can find an industrial 
justification in the market of monitoring equipment. 


1.3 Detection of the Displacement Vector 


The interferometric radar detects only the displacement component along the 
view direction. The real displacement (of slope, building, bridge, tower, ...) 
is retrieved with some hypothesis about the motion. For example, in a slope, 
the displacement can be supposed along the gradient, in the lower deck of 
a bridge along the vertical axis and in a tower along the horizontal. These 
suppositions can be reasonable in most cases, but there are important excep- 
tions. As Dei et al. in 2013 [79] demonstrated for bridges, some targets, that 
are not at the center of the span, have a component of horizontal displace- 
ment during static monitoring of bridges. If this horizontal displacement is 
not considered, the results are absurd (a deck that rise when it is loaded). 
Figure 1.11 shows how this effect is in the relationship with the thickness of 
the deck and the position of its neutral axes. 


Figure 1.11: Effective displacement of a point P at an edge and close to the 
central span during the static test of a bridge. 


The case of slope monitoring can be more complex. For example, for 
detecting the real displacement pattern of Lamosano village in Italian Alps, 
built on a landslide, it has been necessary to combine data from both in- 
terferometric radar and lidar [42]. Unfortunately, this complex behaviour of 
landslides is rather common. 

In order to obtain more than one component of the displacement vector, 
in 2014, Severin et al. [80] used two GBSAR systems in different positions 
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that monitor the same slope. In 2015, Zeng et al. [81] proposed the same 
principle using two MIMO. Two radar systems operating independently can 
be a solution, but surely is expensive and not practical. 

Other solutions are effective only in short range; Hu et al. [82] proposed 
to divide the scan in two sub-apertures for obtaining two points of view 
with the drawback to reduce the angular resolution; Pieraccini et al. [83] 
proposed to use the rotation of the antenna in a vertical plane for obtaining 
two or three components. They tested the technique during the static test 
of a bridge in Florence, Italy [84]. 


1.4 Three-dimensional Imaging 


The idea of obtaining three-dimensional radar images using a GBSAR started 
with this radar itself. Already in 2000 Lopez-Sanchez [85] described a SAR 
migration algorithm for focusing three-dimensional images in an anechoic 
chamber. In the same year, Tarchi et al. [86] developed a GBSAR that 
obtains a 3D image of the reproduction (in scale 1:2) of an ancient facade. 

Generally speaking, complex targets like buildings, tunnels, and construc- 
tion works would require 3D imaging capability, but there are two practical 
problems. The first is the long acquisition time. With reference to Fig- 
ure 1.12, in order to image in the three-dimensional space using a liner SAR, 
the radar head must scan step-by-step and line-by-line a planar surface. 
As an example, a fast GBSAR able to scan a line of 100 points in 5 min- 
utes [87], will scan a plane of 100 lines in more than 8 hours. It means that 
the radar is able to detect only very slow movements. The second problem 
is the size and bulk of the mechanical system, which makes its deployment 
hard. 

These are the main reasons why GBSARs with 3D imaging capability 
are not very popular. 

The use of higher frequency (e.g., W-band) mitigates the problem of the 
bulky equipment, and partially even the problem of acquisition time, as the 
movements are smaller and, therefore, faster. The radar HYDRA [88] is a 
step in this direction. 

A more radical solution is a two-dimensional MIMO [89]. These are 
exciting developments, but they are again prototypes in a laboratory and 
not suitable to be used in the field. 


1.5 Objective 


More than 415 scientific papers on GBRI have been published from 1997. 
The linear SAR is the most common system, followed by the no-DOA sys- 
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Figure 1.12: GBSAR with 3D imaging capability. 


tems. The majority of radar systems work in Ku-band. More than 40% of 
case studies are relative to landslides followed by bridges (30%) and buildings 
(20%). 


There are several open issue in the GBRI systems that can involve their 
application, e.g. wind turbines, or their technological improvement. The 
main technological development directions are forward faster radar (i.e. from 
SF modulation to FM modulation), radar with wider view angle (using Arc- 
SAR configurations [22], [23], [24]), MIMO radar, radar with capability to 
detect the vector of displacement and not only a single component, and small 
size radars operating at higher band (in Ka or W band). 


Currently, FM modulation and ArcSAR appear already mature for the 
market of monitoring equipment. Some manufacturers have products with 
these characteristics. MIMO GBSAR is at a good level of industrialization 
and probably, in a short time, we will see its diffusion. The detection of the 
displacement vector is again at the level of prototypes, but its usefulness is 
evident, and it will probably soon be established as a standard technique. 
Finally, interferometric radar in Ka or W band are again laboratory proto- 
type, but probability they will have a role in the future, as well as radars 
with 3D imaging capability. 
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The aim of this work is to present some solutions founded during the 
doctoral course of the author. In the next chapters a bistatic system for 
retrieving the displacement vector, a Compressive Sensing MIMO system 


and some 3-D radar systems are presented. 
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Basic principles of Ground-Based Radar Interferometer 


In this chapter a theoretical background of GBSAR are introduced. The 
main definitions (i.e. range resolution, azimuth resolution) and concepts 
(i.e. focussing algorithm) of ground-based/terrestrial radar are summarize. 
In the next chapter, when it is needed, the concepts are extended or resumed. 


2.1 Target detection 


Radar is an electromagnetic devices able to locate reflective objects called 
targets. The basic principle of radar is illustrated in Figure 2.1. A transmit- 
ter generates a signal (e.g. a short pulse or a sine function) that is radiated 
by the TX antenna. A portion of the transmitted signal is backscattered by 
the target in many directions. The backscattered signal directed back to- 
wards the radar (named echo) is collected by the RX antenna which delivery 
it to a receiver. Therefore it is processed to detect the presence of the target 
and its location. 

General speaking, the distance, or range, to a target can be calculated 
by measuring the time that signal takes to travel to target and return back 
to the radar, according to the follow equation: 


c AT 


R= 5 


(2.1) 
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Figure 2.1: Basic principle of radar. 


Where R is the range to the target, c is the speed of light and AT is the 
travel time. 

Radars can operate by using several kinds of signal: short pulse, Stepped- 
Frequency, frequency modulation. Usually, Stepped-Frequency (SF) wave- 
form (see Fig. 2.2) is used in terrestrial interferometer. Indeed, the first 
terrestrial interferometric radar used a Vector Network Analyzer (VNA) as 
transceiver [5] that operated a Stepped-Frequency (SF) waveform. 

Most of the commercial equipments use SF modulation, even if, recently 
the radar interferometers of IDS geo-radar [90] and meta-sensing [91] are 
both Frequency-Modulated (FM) systems. 


In this thesis most of the radar systems operates a SF signal (after [24]). 


Figure 2.2 shows an example of SF waveform. The SF radar sends a 
sweep of frequencies from the starting frequency fı to the end frequency 
fn; = fı + B, where B is the bandwidth, with Ny steps. 


Inv 


Figure 2.2: Step Frequency modulation. 


The in-phase (/;) and quadrature (Q;) component of scattered signal is 
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detected for each j-frequency 
Ej = 1; + iQ; (2.2) 


where E; is echo received by radar and è is the imaginary unit. 


The radar image can be obtained by calculating the inverse Fourier trans- 
form (IFFT) of the echo [24]: 


I =IFFT(E, Ny: F) (2.3) 


where I is the obtained radar image and F is the padding factor. 
In order to reduce the side lobes and improve the quality of the image a 
frequency windowing can be applied to the echo before calculating the IFFT 


EY = Ej Wj (2.4) 


where W is the windowing function and EW is the j-component of echo after 
window. 

There are several windowing functions, Figure 2.3 shows an example of 
some possible of them [92]. The windowing decrease both the side lobe level 
and resolution. 
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Figure 2.3: Example of windowing functions. 
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In this thesis a Kaiser window has been used: 


» (oi @ a) 


1- 
Io(8 


w(n) = (2.5) 


where 0 < n < N, N is the number of samples of the window (equal to the 
number of frequencies), Io is the zeroth-order modified Bessel function of the 
first kind and £ is a non-negative real number that determines the shape of 
the window. The ( used in this thesis is usually equal to 5. 


Range resolution and unambiguous range 


The transmitted signal utilized a frequency bandwidth in order to separate 
and distinguish the targets. Indeed, single targets can be discriminated in 
range if their separation is bigger than [93]: 


Tes — 35 2.6 
R 5B (2.6) 


This equation represents the range resolution (Rres) of a radar. 

The unambiguous range of a radar is the maximum range at which a 
target can be located so as to guarantee that the reflected signal/pulse from 
that target corresponds to the most recent transmitted pulse. In other word, 
if a target is farthest than the unambiguous range, it appear to be close to 
radar system as shown in Figure 2.4. 
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Figure 2.4: Unambiguous range. 


In reference with Figure 2.4 the unambiguous range, Rmaz, can be found 
as follow: 


N; (2.7) 
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Azimuth resolution and angular ambiguity 


While the range can be detected by using the bandwidth of the transmitted 
signal, the direction of arrival can be detected by exploiting the movement 
of the antennas. 

The azimuth resolution is the minimum angle separation in order to 
resolve different target located at the same range (see Figure 2.5(a)). If 
the angular separation of two targets in the same range is smaller than the 
azimuth resolution (full width half maximum), their coherent summation 
results to be one single peak as shown in Figure 2.5(b). 


arget 2 
i Sum 
Resolution 


(a) Resolved target (b) Non-resolved target 


Target 1 


Figure 2.5: Azimuth resolution. 


The most advanced radar systems are able to detect the direction of 
arrival through synthetic aperture (SAR: Synthetic Aperture Radar). In 
Figure 2.6 is shown the basic principle of a linear GBSAR, but all the followed 
conclusion can be generalized for any GBSAR system. 

With reference to Figure 2.6(a), the signal received by the radar in posi- 
tion x=0, at central frequency fo, from a target located in Ro is 


E = ei foPo (2.8) 
while, the signal received in a generic position x is 
E(x) = ei E fo(Ro—# sin 00) (2.9) 


A rather easy method for obtaining a spatial image of a GBSAR is to 
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Figure 2.6: Basic principle of a linear GBSAR. 


calculate the Fourier transform of the signal in spatial domain. 


f(0) zS E(x)et Ersin b) dy — 


L J-r/2 
2.10 
e tS fo Ro fe An P Ag . ( ) 
= Ji et E (fox sin 00) oj F (e sin 0) da 
L - 1/2 


where f(0) is the radar image. 
By solving the integral in (2.10), f(0) is 
._ | QL fo sos i 
sin [2 (sino — sin 00) 
0) = e` iE foRo c 2.11 
f(O) 2rLfo (sin 9 — sin 0o) CI 


The Equation (2.11) represents a sinc. Conventionally the full width half 
maximum (resolution) of a sinc is given 


2r L fo 


(sin — sin 0o) = 7 (2.12) 


The angular resolution can be found by inserting the Taylor series ex- 
pansion of sin 6|g=9, in (2.12) 


4TLfo 


cos 0oA0 ~ 1 (2.13) 
where the second order of the series was used. 
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Therefore the angular resolution A0 is 


À 
Ad S ——_ 2.14 
2Lcos0o ( ) 
with A wavelength. The definition in (2.14) can be generalized for any dif- 


ferent Synthetic Aperture system. 


Angular ambiguity 


With reference to the unambiguous range, it is possible to define an angular 
ambiguity for the azimuth direction. 

Generally speaking, to avoid the angular ambiguity, the step between 
two position in the GBSAR systems has to be smaller then a quarter of 
wavelength. 

As shown in Figure 2.6(b) the phase difference between two radar posi- 
tions, with step Ax, of a target in a generic angular position 0 can be write 
as Az sin(6). 

The phase difference between to position has to be 


2x Axsin(0) = 2 (2.15) 


in order to focus the target in the same position as shown in (2.9). 
If Ax = \/2 the ambiguity, 0u, is in +7/2 (Figure 2.7(b)). Instead, if 
Ax < A/4 all possible angular ambiguity is avoided. Figure 2.7 shows an 


example of angular ambiguity. 

In real application the constrain given by (2.15) can be relaxed. Indeed, 
in real application the antennas are usually directional and it is enough that 
Ou is considerably larger than the angular aperture of the antenna. 


2.2 Focusing algorithm 
The measured echo signal in (2.9) can be write as follow: 
Ejk = ljk + iQj,k (2.16) 


where E; x is the echo of j-frequency taken from k-position of the radar. 
With reference to Figure 2.8, the basic formula for focusing on a generic 
n-point (x, y, z) of the space surrounding the radar is [24]: 


in Hype e Pei (2.17) 
j,k 
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Figure 2.7: Example of angular ambiguity. 


Equation (2.17) is rather computational heavy, so the key point of any 
focusing algorithm is to reduce its computational cost without a significant 
loss of generality and accuracy. A simple focusing algorithm has been pro- 
posed by Pieraccini et al in 2017 [24] 


I(x,y,2) = $ Ah @-DAREPT (Ey, Ny- Fete FA-@ YAR) (2.18) 
k 


where I(x, y, 2) is the radar image in point (x,y,z), c is the speed of light, fi 
is the fist frequency used in SF signal, B is the bandwidth, N} is the number 
of frequency, F is the padding factor, AR is define as 

T N fT lec 


AR= 2.1 
R= NF 2B n. 


and p is an index defined by the k-position of antenna n-position on the grid: 


p = integer (Fe) (2.20) 
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Figure 2.8: Focusing grid. 


This focusing algorithm can be generalized for any kinds of synthetic 
aperture and signal used. 

To improve the image quality a window can be applied also in cross-range 
direction. 


2.3 Interferometry 


The interferometric radar uses the phase information of the received signal 
to measure small displacement of the targets. Figure 2.9 shows the working 
principle of this technique, detailed description can be founded in the [8], 
[94], [95]. 


Figure 2.9: Working principle of ground based radar interferometry. 
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With reference to Figure 2.9 the phase difference between two acquisitions 
is given by: 
Ay = Apdisp + A Patm + A Pnoise (2.21) 


where Agdisp is related to a possible displacement of the target, Ayatm 
is related to the atmospheric conditions (humidity, temperature, pressure), 
and AYnoise is a phase noise contribution that can never be completely 
eliminated in experimental measurements. 

The atmospheric phase contribution (usually named “atmospheric phase- 
screen” [96]) can be neglected for short-range short-term applications (in- 
dicatively when the distance radar-target is less than 100 - 200 m and the 
time between the acquisitions is shorter than 1 - 2 hours). In all other 
cases, the atmospheric phase-screen has to be removed with suitable pro- 
cessing techniques. The proposed algorithms for atmospheric correction are 
many [97], [98], [99]. Their working principle is always based on the identi- 
fication of stable targets in the radar image. Using these stable targets the 
phase is compensated using a suitable (always approximated) atmospheric 
model. This phase-screen removal is a very critical point and limits the final 
accuracy in long term monitoring. 

When the phase noise contribution is negligible and the atmospheric 
phase-screen is removed (or negligible), the possible displacement (AR) of 
the target between two acquisitions is detected as phase difference (Ay) 
through the following basic relationship [100]: 


A 
_ Av) 


AR= iz 


(2.22) 
It is important to note that an interferometric radar is able to detect only 
the component of displacement in range direction. 

Generally speaking an interferometric radar can be used both for dynamic 
and static test. 

Dynamic Monitoring: in this case the main interest is to be able to 
retrieve the frequency spectrum of displacement. Nowadays, the dynamic 
monitoring is suitable for slander structure because of the high sample fre- 
quency (more than 100 Hz) does not allow to move the radar head to achieve 
the azimuth resolution. 

Static Monitoring: the movement in this case is very slow (quasi- 
static) and the main interest is to detect precisely the displacement. While 
the frequency spectrum measurement does not present critical issue, the 
measure of absolute displacement is affected by systematic and determinis- 
tic error (phase noise, phase-screen, movement of structures, ...) that can 
be kept small but cannot be entirely eliminated [101]. In this case the mea- 
surement time is long enough to allow the movement of the radar in order 
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to achieve azimuth resolution. The static monitoring can be used for infras- 
tructure (bridge and building) or for natural environment such as landslides. 


In the next Chapters some prototypes of advanced interferometric radar 
are presented. Unfortunately, most of the measurements presented are re- 
lated to static monitoring. This is due to the transceiver used: an old model 
of VNA (HP8720A). This VNA had a time-constraint for the data transmis- 
sion with computer. Indeed , some of the radars presented could be used 
also for dynamic monitoring by using a more recent equipment. 
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Radar systems for retrieving the displacement vector 


The GBRI systems are able to provide maps of displacement along the range 
direction. Nevertheless, Severin et al. [80] demonstrated the great advantage 
of the vector displacement in slope monitoring, and Dei et al. [79] experimen- 
tally showed how the detection of only one component can give paradoxical 
results in structural monitoring (a bridge deck that is apparently raised when 
loaded). 


In 2014, Crosetto et al. [102] proposed a non-interferometric procedure 
for detecting the transversal component of the displacement. In 2017, Hu 
et al. [82] proposed to process two half scans separately for obtaining two 
images with a baseline. Both methods are limited by the scan length: the 
larger is the scan the better is the accuracy of the traversal displacement. 


In this chapter a GBSAR operating in monostatic and bistatic modality 
is presented [103], [104]. This system made use of a transponder as addi- 
tional channel for performing the measurement. In this way, the same radar 
equipment is able to acquire two images taken from different points of view: 
so it is able to acquire two different components of the displacement of the 
targets. Contrary to the works [102], [82], the accuracy in the detection of 
the transversal displacement does not depend on the scan length (that in 
a GBSAR, installation can hardly overcome 2 - 3 m) but on the distance 
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between radar and transponder that can be tens or even hundreds meters. 
In the last part of the chapter some other applications of this principle 


are presented [105], [106]. 


3.1 Working principle of bistatic 

The working principle of the bistatic GBSAR is shown in Figure 3.1. A linear 
GBSAR acquires an image of its own field of view. Using a third antenna, 
it acquires a second image of the same scenario exploiting the bounce of 
the signal through a transponder. The transponder consists of a couple of 
antennas and an amplifier system. 

The bistatic SAR image can be focused through the backprojection al- 
gorithm explained in Chapter 2. With reference to Figure 3.1, the back 
projection is implemented by compensating the phase history of any path 
radar-target-transponder-radar [104]. 

The monostatic and bistatic images are focused in the same grid, but 
generally this does not assure that they are coregistere [107], so they could 
need the application of a specific processing technique for associating the 
corresponding pixels in the two images. Nevertheless for targets that are 

smaller than resolution, and in an environment without heavy clutter, we 
can assume that the two images are coregistered. The experimental setup 


has been arranged according to this assumption. 


Target 


n RX 
it E f 
it 
re 
d TX Transponder 


Figure 3.1: Working principle of bistatic GBSAR. 


As known, in the monostatic configuration, when a target in the position 
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(xo, yo) is displaced of a vector 5, the electromagnetic path is changed of 
Aly = 2(8- ü) (3.1) 


were Ñ is the unit vector from the radar to the target position (see Figure 
3.2). 


Figure 3.2: Detected displacements. 


In accordance with (2.22) of the fundamental interferometric relation 
between phase (Ag) and displacement (Al) we can define: 


Lu A 
Ası = 8- Ñ = gAn (3.2) 


with Ay, differential phase measured by monostatic configuration. 
In the bistatic configuration, when a target in the position (xo; yo) is 
displaced of a vector 5, the electromagnetic path is changed of 


Al = f- ti, + 8: iy = 8. Hi + | =28- Boos (5) (3.3) 


where Ñz is the unit vector shown in Figure 3.2, b is the bisector unit vector, 
and £ is the bistatic angle. 

Therefore the displacement measured by the bistatic configuration 
À 1 
4r cos( 6/2) 
with Ag»ə differential phase in bistatic configuration. The detected displace- 
ments Ası and Aso are in linear relationship with the two orthogonal com- 
ponents Ax and Ay of the effective displacement 


(1) (3) es 
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with 
a oe o 

and 
cos Vp ™ = e : m= 1,2 (3.7) 
cos By" = i m=1,2 (3.8) 


where è and j are the orthogonal basis of the Cartesian plane. 


Therefore 
Aa i As, 
= M7! f 
(S) (An) na 


Equation (3.9) allows to calculate Ax and Ay from the measured As; 
and Aso. 
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Figure 3.3: Measurement geometry. 
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In order to assess the measurement error, we can consider the case shown 
in Figure 3.3. The target is in front of the radar (along the y-direction) and 
the transponder is on a side. For this geometry, the matrix M results 


0 1 
As eB (£) cos 2) (30) 


Let dy, be the phase uncertainty in monostatic configuration and 692 be 
the phase uncertainty in bistatic configuration, by inverting the matrix M 
and using (3.2) and (3.4), the displacement uncertainty along y (dy) results 


= 


ô 
Y 4T 


dpi (3.11) 


while the displacement uncertainty along x (dx) results 


COS 2 An Yi 


This equation expresses how the uncertainty along x-axis depends on 
bistatic angle 6. It is worth to note that dy, and dy can be different. 
Indeed, dy2 depends on bistatic radar cross section of the target (that is 


dr = 


2 | A 
Ta i (512) 


different from monostatic radar cross section) and even on the gain, noise 
figure, and phase stability of the transponder. 


3.2 Range and angular resolution of bistatic GB-SAR 


For estimating the range resolution we can refer to Figure 3.4. Considering 
the path C-A-T and the path C-B-T, they are separated in time if their differ- 
ence is larger then the time resolution ¢/B, with c speed of light and B band- 
width. 
Therefore 
AR+ ARcos 8 = 5 (3.13) 


that is 


yB 


PAs 2 - cos? (5) 


(3.14) 


where ĝ is again the bistatic angle. 

It is worth to note as it differs from the well-known range resolution for 
a monostatic radar only for its dependence on the angle 8 and it coincides 
with the later when 8 = 0. 

For estimating the angular resolution we have to take into account that 
the locus of points at the same time-delay are the points of the ellipse that has 
the scan center (0,0) and the transponder (a7, yr) as foci (see Figure 3.5). 
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Figure 3.4: Estimation of range resolution. 


Therefore it is convenient to define the angular resolution Age with respect 
to the point (ac, yc) along the bisector of angle 6 crossing the straight line 
joining the scan center (0,0) and the transponder (x7,yr). If we shift the 
origin of axes in the point (xc, yc) the angular resolution will result well- 
defined. 

By geometric consideration: 


to = Recos(y) yo = Rc cos(y) (3.15) 


where Rc is the distance of (xc, yc) from the origin of axes. It can be 
obtained as 
sin (3) 


Ro = |— ~] Ro (3.16) 
cos Q +9- 3) 
with Ro in Figure 3.6. 
Figure 3.6 sketches how to estimate the quantitative value of the angular 
resolution Agc. 
It results with elementary geometrical considerations 


ADR 


= A .1 
by using Ad’ from (2.14), Adc results 
Adc Ro a (3.18) 


~ Ra cos(8/2) 2L cos 3 
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Figure 3.5: Constant time-delay ellipse. 


3.3 Proof of principle 
Equipment 


Figure 3.7 shows a sketch of radar prototype that has been assembled for 
testing the working principle above described. The radar front end has been 
fixed on a mechanical axis controlled by a step motor. The maximum length 
of the axis was L = 1.80m. 

A vector network analyzer (VNA) HP8720D operated as transceiver pro- 
viding a continuous wave stepped frequency signal in the X-band with cen- 
tral frequency fe = 10 GHz, bandwidth B = 160 MHz, and number of tones 
N; = 401 + 801. The radar operates step-by-step both in positions and 
frequencies. 

A personal computer controlled the axis position and the data commu- 
nications with VNA. 

In reference to Figure 3.8, two RF cables linked the VNA to the front 
end. The VNA output power is 0 dBm, the one-way cable loss is -5 dB, 
and the TX amplifier gains 10 dB. Two single-pole double-throw (SPDT) 
switches provide a direct path (through a 40 dB attenuator) between the 
transmitter and the receiver in order to perform calibrated measurements. 
The aim of the calibration is to avoid the movements of the cables can affect 
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Figure 3.6: Geometry for the calculus of the angular resolution. 


the phase measurement. An other SPDT switcher select the antenna that 
receive the echo from environment (monostatic channel) and the antenna for 
the signal trough the transponder (bistatic channel). The gain of the RX 
amplifier is 20 dB. The integration time of each tone was 1 ms. 

The transponder consisted of two horns with 18 dB gain and a wideband 
amplifier (nominal gain: 50 dB, band: 6-18 GHz, and noise figure: 5 dB) on 
a tripod (see Figure 3.9). 

The two antennas could be oriented and they were separated by a length 
of 1.44 m for minimizing the coupling and the possible oscillation of the 
transponder. Further details of transponder, i.e. its oscillation, will be given 
in the next section. 


Experimental test 


For testing the capability to obtain bistatic radar images, the GBSAR and 
the transponder were installed in a garden of the University of Florence 
(Figure 3.10). At about 35 m in front of the radar there was a metallic 
fence, on the right a line of trees and a metallic signpost. The parameters 
used were: L = 1.80 m, N, = 240, fe = 10 GHz, B =160 MHz, N; = 801. 

The transponder has been installed on the x-axis at 10 m on the left of 
the radar. 
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Figure 3.8: Block scheme of radar. 


Figure 3.11 shows the focused bistatic image. The fence is well visible at 
about 35 m and on the right we can recognize the signal from the metallic 
singpost. 


Figure 3.12 shows the plot of the amplitude in function of de (the angle 
with respect to (xc,yc)) at the range of the metallic signpost that gives 
the maximum signal in the bistatic image. The plot relative to the target 
has been compared with the simulation of a point-target. The two plots 


were normalized at their maximum. The agreement between measurement 
and simulation is excellent and it confirms the bistatic radar imaging theory 


previously described. 
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Figure 3.9: Transponder. 


Radar 


Transponder 


Figure 3.10: Aerial view of the University’s garden. 


Measuring two components of displacement 


The capability to measure two components of displacement has been tested in 
a controlled environment by making use of the same setup described before. 
The test site is an open field of University in Sesto Fiorentino, Florence, 
Italy. 

A special target has been assembled aiming to test the capability of the 
radar to detect two components of the displacement. The target is a metallic 
pole provided with a small adjustable corner reflector 10 x 10 cm? and three 
grids disposed orthogonally, as shown in Figure 3.13. The aim of the grids 
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Figure 3.11: Bistatic radar image [104]. 
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Figure 3.12: Measured angular plot of the metallic signpost compared with 
the simulation of a point-target [104]. 


is to provide a significant bistatic radar cross section. 

The pole has been mounted on a linear micrometric stage platform with 
two orthogonal controlled movements on the xy plane with 0.1 mm nominal 
accuracy and 25 mm of maximum displacement. Some eccosorb has been 
used to cover the mechanical device. The micrometric stage was fastened to 
a wood table and a dead load kept the wood table fixed. 

The target and the transponder were positioned, as shown in Figure 3.3. 
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Figure 3.13: Picture of the target [103]. 


The distance target radar (yo) was 104.45 m, and the coordinates of the 
transponder were xr = 39.53 m and yr = 45.80 m. 

Figure 3.14 shows the obtained monostatic image of the scenario where 
the radar has been installed, compared with an aerial picture. The more 
evident features are the large building at 300 m in front of the radar, the 
shack at about 200 m, and the corner of the building at about 130 m on the 
right. The target appears as a single point. 


(ap) 1omog 


Figure 3.14: Monostatic radar image [103]. 


The obtained bistatic radar image of the same scenario is shown in Fig- 
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ure 3.15(a). The target is very recognizable as well as a feature (an artefact) 
that apparently does not correspond to any physical target. 

This artefact is due to the clutter behind the transponder that can give a 
monostatic signal even in bistatic configuration. This has been confirmed by 
focusing in bistatic modality an acquisition carried out when the transponder 
was switched OFF (Figure 3.15(b)). 
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(a) Bistatic radar image [103]. (b) Bistatic radar image when the 
transponder was switched OFF [103]. 
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(c) Radar image obtaining, focusing with (d) Bistatic radar image obtaining as 
the monostatic algorithm, the bistatic subtraction of the images in Figure 
data acquired when the transponder was 3.15(a) and 3.15(b) [103]. 


switched OFF [103]. 
Figure 3.15: Bistatic radar image and its artefacts [103]. 


By focusing the same acquisition in monostatic modality (performed 
when the transponder was switched OFF), it is evident that the artefact 
in bistatic image is due to the building behind the transponder on the right- 
hand side of the field of view (see Figure 3.15(c)). 

This suggests a possible procedure for rejecting the artefact: to perform 
a bistatic acquisition when the transponder is switched ON and a second 
bistatic acquisition when the transponder is switched OFF. In the subtrac- 
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tion of the two images, the artefact results just a weak residual, as shown in 
Figure 3.15(d). 

In order to evaluate the capability of the radar to detect two components 
of the displacement, the micrometric stage platform carried out displace- 
ments along x (toward the positive direction of the x-axis) at steps of 5 mm, 
from 0 to 25 mm. The Ax and Ay measured by the radar are shown in 
Figure 3.16. 

The obtained mean absolute deviation of Az around the nominal value 
was 1.74 mm, while the deviation of Ay was 0.60 mm. By considering each 
single step as an independent measurement, the mean value along x was 
5.59 mm (to be compared to the nominal value 5 mm) with a standard 
deviation of 0.84 mm, while the mean value along y was -0.17 mm (to be 
compared to the nominal value 0 mm) with a standard deviation of 0.43 mm. 
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Figure 3.16: Displacements along x. Circles: measured displacements along 
x. Triangles: measured displacements along y. [103] 


In a second measurement session, the same micrometric stage platform 
carried out displacements along y (toward the radar). In this case the dis- 
tance target-radar (yo) was 75.39 m, and the coordinates of the transponder 
were £r = 35.10 m and yr = 37.53 m. The Ax and Ay measured by the 
radar are shown in Figure 3.17. 

The obtained mean absolute deviation of Ax around the nominal value 
was 2.16 mm, while the deviation of Ay was 1.67 mm. By considering each 
single step as an independent measurement, the mean value along x was -0.41 
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mm (to be compared to the nominal value 0 mm) with a standard deviation 
of 3.02 mm, while the mean value along y was 5.33 mm (to be compared to 
the nominal value 5 mm) with a standard deviation of 0.64 mm. 


Nominal displacement (mm) 
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Figure 3.17: Displacements along y. Circles: measured displacements along 
x. Triangles: measured displacements along y. [103] 


The capability of a monostatic/bistatic GBSAR to acquire the displace- 
ment vector was demonstrated with these two measurements section . 

Some critical issues can be found in this basic approach: gain of transpon- 
der, “gosth targets”, long-term phase stability of transponder. 

In the next Section possible solutions of these issues will be presented. 
Anyway, for some of them a simple solution can be found by using the before 
setup. 

For avoiding the presence of “ghost targets” in the bistatic image, as a 
practical rule, the transponder should not be installed in a position, where 
there is heavy clutter in the background. Nevertheless, an effective coun- 
termeasure for rejecting/mitigating the “ghost targets” was to subtract the 
image acquired when the transponder was OFF from the image acquired 
when the transponder was ON. Anyway, the scenario can have small changes 
between the two acquisitions, so a weak residual could remain in the bistatic 
image. 

A possible phase-changing of the transponder (phase of amplifier, small 
uncontrolled movement of the mechanical support of the transponder) can 
be treated as a change in the atmospheric phase screen and so can be com- 
pensated. 
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3.4 The transponder 


As before mentioned a basic transponder consists in a couple of antennas 
and an amplifier. The gain of the amplifier has to be enghout to compensate 
the path-loss and the noise. 

In reference with Figure 3.18, by applying the radar equation [93] in 
monostatic and bistatic configurations, we obtain the following equations, 


respectively: 
p(mono) = P,.G 1 (mono) È G `? (3 19) 
cala 4nR? dn 
i de i ae 1 A? 
p(rist) — PGi (brst) Ci — AG, Gs (3.20) 


mR, mR Sin a RE, a 


where Py is the transmitted power; G1, Go, G3, G4, Gs are the gains 
of the five antennas; 1, is the distance between the radar head at the 
k-position along the mechanical guide and the target; Ri is the distance 
between the center of the mechanical guide and the target; Rə is the distance 
between the target and the transponder; A is the gain of the amplifier of the 
transponder; Rs, is the distance between the radar head at the k-position 
along the mechanical guide and the transponder; R3 is the distance between 
the center of the mechanical guide and the transponder; 0(mono) is the radar 
cross section in monostatic configuration; 0(yst) is the radar cross section in 
bistatic configuration; and A is the wavelength. 
In practical deployments, the following assumptions are reasonable: 


1. All the antennas have the same gain G, = G2 = G3 = G4 = Gs 
G. 


2. The length (Lo) of mechanical guide is much smaller than the radar- 
target and target-transponder paths > Rx = Rı and R3,x = R3. 


3. Although R1 and R3 are not equal, they can be reasonably similar > 
Ri = R3. 


4. (mono) ~ F(bist)- 


The last assumption is rather hard, but it can be done when the aim of 
the calculation is just a rough estimate of the received power. With these 
assumptions, the (3.19) and (3.20) become: 


(3.21) 


poro SA GA 2 
piena) E 4r R3 
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Figure 3.18: Monostatic/bistatic GBSAR with transponder. 


Hence, in order to compensate the loss due to the path transponder-radar, 
the gain of the transponder should be as follows: 


Ar R3 $ 

A= ( GI ) (3.22) 

This gain value compensates the loss; however, with this value of gain, 
the noise figure of the transponder is added to the noise figure of the first 
amplifier of the radar head. In order to improve the noise figure, the gain 
has to be increased. Indeed, as it is well known, the noise figure of the whole 
system is equal to the only noise figure of the transponder amplifier (i.e., the 
first stage of the receiver chain) if the gain of this is much greater than the 


loss of the path. In practice, a reasonable rule of thumb for the gain of the 
transponder could be the following: 


Lol (3.23) 


A>l 
< o( GX 


For obtaining good displacement measurements in both directions with 
the bistatic configuration, L should be in the order of magnitude of the 
distance between radar and target (see (3.12)). 

The operative distance of a GBSAR can reach up to 5 km. However, a 
value between 50 m and 500 m is reasonable in most applications. Therefore, 
by considering an antenna gain (G) of 15 dB and a wavelength of 3 cm (10 
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GHz), we obtain that the transponder gain (A) could range between 66 dB 
and 86 dB, with L ranging between 50 m and 500 m. 

However, a major problem of any transponder is its possible oscillation. 
The isolation between the two antennas has to be higher than the gain in 
order to avoid the oscillation. 

A way to increase the isolation between two linearly polarized antennas 
is to rotate one of them by 90°. For real antennas, the cross-polarization can 
provide an isolation of 20-30 dB in addition to the co-polarized isolation (see 
for example the reference in [108]). In the specific case of the transponder 
described above, the cross-polarization of the antenna pointing at the radar 
does not affect the bistatic operation if the second receiving antenna of the 
radar head is cross-polarized with respect to the transmitting antenna (see 
Figure 3.19). 


O, Target 


RX 


< ; Transponder 


HH TX 


Figure 3.19: Monostatic/bistatic GBSAR with cross-pol transponder. 


A further advantage of the cross-polarization is that monostatic and 
bistatic signals have different polarizations, so possible artefacts due to the 
mixing of bistatic are greatly mitigated. 

This simple solution is not enough in practical application when the 
distance radar-transponder has to be higher than 100 m. 


Cross-pol Transponder with frequency shifter 


Another effective way to avoid the oscillation of a high gain transponder is 
to shift the output frequency using a suitable frequency shifter. Figure 3.20 
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shows the scheme of a bistatic radar that uses a transponder with frequency 
shifter. 


Figure 3.20: Monostatic/bistatic GBSAR with frequency shifter transpon- 
der. 


The radar operates at central frequency fo. The transponder has its own 
local oscillator (LO). A mixer and a high-pass filter operates as up-converter 
and shifts the central frequency, so the frequency of the signal transmitted 
from the transponder to the radar head is fo + fro. In the radar head, 
a down-converter (a mixer and a low-pass filter) brings the signal back to 
frequency fo. The two LOs have to be coherent and synchronized, so a 
reference signal at low frequency frer is shared with a cable. 

As the cable could move during the movement of the radar head along 
its mechanical guide, a calibration path is necessary to compensate possible 
phase shift of the reference signal. Moreover, the calibration provides long- 
term phase stability and compensates possible phase changes that can occur 
when the equipment is switched off and on again. Therefore, the radar head 
has a second transmit channel (selected by a Radio Frequency switch) that 
provides a calibration signal directly transmitted to the transponder. The 
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transponder has two receiving antennas (one pointed at the target and the 
other one pointed at the radar). A RF switch, controlled by a Wi-Fi device, 
selects between the two antennas (measurement and calibration). 


Test Bench of the System 


Figure 3.21 shows the implementation of a up-converter realized with a mixer 
(SPECWAVEMM96P-L15), a LO at 2.0 GHz (HP 8672A), and a high-pass 
RF filter (FILTRONIC SB 029). Two band-pass filters (MINICIRCUITS 
VHF-3800+ and FBP-5.8-11G tuned in the band 5.8 GHz - 10.1 GHz) are 
used for limiting the input band. 


Signal Generator HP 8672A 


pa SPECWAVE i it 
MINICIRCUIT FBP-5.8-11G MMSGR:LIE i FILTRONIC | 


10 GHz- | VHF-3800+ Tuned 5.8-10.1G po SEU ne 19 GH 


Figure 3.21: Up-converter operating at 10 GHz with 2 GHz frequency shifter. 


The isolation of the transponder with the up-converter was measured 
with the test bench depicted in Figure 3.22. 

Figure 3.23 shows the amplitude plot in frequency of the signal before 
the up-converter (“direct path” in Figure 3.23) and after. The minimum 
isolation resulted in 71.90 dB, which was enough to prevent the oscillation 
of an amplifier with gain up 50-60 dB. 

The up-converter requires a down-converter at the radar head. In order 
to test the capability of the up and down-converters to correctly detect the 
phase of the radar signal, the test bench shown in Figure 3.24 was imple- 
mented. In order to simulate the measured differential phase of a target, a 
manual phase shifter (ARRA 9428B) was inserted between the Vector Net- 
work Analyzer (VNA) and the up-converter. The two signal generators at 
2.0 GHz were locked up with a 50 m coaxial cable that shared the 10 MHz 
reference signal. 

Figure 3.25(a) shows the phase measured when the signal passes through 
the up and down-converters and the phase measured directly (bypassing the 
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Figure 3.22: Test bench for measuring the isolation of a up-converter oper- 
ating at 10 GHz with 2 GHz frequency shift. 
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Figure 3.23: Plot in frequency. Red represents the amplitude of the signal 


through the chain with up and down-converters; blue represents the ampli- 
tude of the signal that bypass the chain through the “direct path” [109]. 
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Figure 3.24: Test bench for testing the phase shift detection. 


up and down-converters). The x-axis represents the marker number of the 
manual roll of the phase shifter (ARRA 9428B). 


Figure 3.25(b) shows the plot of the differences. The average difference 
was 2.78 deg. This difference appears to be a systematic error (probably due 
to nonideality of the mixers and filters), which is compensated by the calibra- 
tion procedure and which cannot anyway affect the differential measurement 
of the targets’ displacements in any case. 


phase [deg] 
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(a) Measured phase by rotating the han- 
dle of the manual phase shifter. Red 
represents the phase measured through 
the up and down-converters; blue repre- 
sents the phase measured bypassing the 
up and down-converters. 
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(b) Difference between the phase mea- 
sured when the signal passes through the 
up and down-converters and the phase 
measured directly (bypassing the up and 
down-converters). 


Figure 3.25: Phase phase measured when the signal passes through the up 
and down-converters and the phase measured directly. 
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Field test 


Figure 3.26 shows the aerial picture of the test site where the radar equip- 
ment was installed. The radar was facing a seven-storey building at 130 
m distance. The transponder was at the right side of the radar at 22.2 m 
distance. Figure 3.27 shows a picture taken by the radar position. 


Figure 3.26: Aerial view of the test site [109]. 


The radar operated in the band 9.915 GHz - 10.075 GHz. The number of 
frequencies was 801. The length of the mechanical scan (orthogonal to the 
view direction) was 1.82 m. The number of steps along the scan was 180. 
The transmitted power was 19 dBm. The measured gain of the transponder 
(shown in picture in Figure 3.28) was 91 dB. 

Figure 3.29 shows the obtained monostatic image. The edges of the 
building marked in the pictures in Figures 3.26 and 3.27 are well evident. 
The amplitude scale of the image is in signal to thermal noise ratio (STNR), 
which has been calculate as outlined in Reference [110]: 


2 
Po TmeasAcal 


STNR;; = zi (3.24) 
B 


where i °°) is the radar image focused using the algorithm; P) = 19 dBm is 
the transmitted power; Tmeas is the integration time that is given by the time 
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Figure 3.27: Picture taken from the radar position in front of the building 
in the field of view of the radar [109]. 


Figure 3.28: Picture of the transponder in the field [109]. 


of a single continuous wave (CW) measurement multiplied by the number of 
frequency (N; = 801) and the number of positions (N, = 180); Acar = 40 
dB is the attenuation of the internal calibration path of the radar head; kp 
is the Boltzmann’s constant; T is the absolute temperature (300 K). 

In order to test the phase stability of the radar system, 13 images in 7 
h were acquired. The interferograms has been calculated between any image 
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Figure 3.29: Monostatic radar image [109]. 


and its subsequent image. Figure 3.30 shows one of these interferograms. 
The pixels relative to the building appear rather stable in-phase. 
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Figure 3.30: Interferogram between two monostatic radar images [109]. 


The obtained bistatic image is shown in Figure 3.31. The image quality 
is comparable with the monostatic image. 

The amplitude scale of the image in Figure 3.31 is STNR calculated 
using (3.24), where Aca; is not longer the attenuation of the internal cali- 
bration path but is the attenuation of the external bistatic calibration path 
estimated as follows: 


2 
A Se GTX-radar GRx-transp»1 A 
cal T att” 


4n R? An (3.25) 
GTX-transp GRX-radar AG A i 
4r R2 An RX—radar 
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Figure 3.31: Bistatic radar image [109]. 


With reference to the block scheme in Figure 3.32, Grx-radar = 14 dB is 
the gain of the transmitting antenna of the radar pointed at the transponder; 
R3 = 22.2 m is the radar-transponder distance; Grx-radar = 14 dB is the 
gain of the receiving antenna of the transponder pointed at the radar; A, 
is the wavelength at 10 GHz; Ag = —30 dB is the attenuation of the 
attenuator in the calibration channel of the transponder; A = 91 dB is the 
gain of the transponder; Grx-transp = 14 dB is the gain of the transmitting 
antenna of the transponder pointed at the radar; Grx_transp = 14 dB is the 
gain of the receiving antenna of the radar pointed at the transponder; Az is 
the wavelength at 12 GHz; and Arx_ radar = —2 dB is the gain of frequency 
shifter of the radar. 

Even in bistatic configuration, 13 images were acquired. Figure 3.33 
shows one of the interferograms obtained with a couple of calibrated bistatic 
images. 

The point marked with “A” was used to quantify the phase stability both 
in monostatic and bistatic interferograms. The “apparent” displacement 
(AR) was calculated as (2.22). 

As the point marked with A is presumably a stable point, any displace- 
ment fluctuation is due to the radar system, atmospheric fluctuation, and 
environmental temperature changes. Figure 3.34 shows the time series of the 
three measured displacements of the point A: monostatic, bistatic without 
calibration, and bistatic with calibration. 

The monostatic displacements had mean value of -0.1238 mm and stan- 
dard deviation of 0.2753 mm. These values are in agreement with other 
measurement campaigns with similar equipment [7], [111], [112]. 

The calibrated bistatic displacements gave mean value of -0.1361 mm 
and standard deviation of 0.3378 mm, in perfect agreement with monostatic 
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Figure 3.32: Bistatic block scheme. 
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Figure 3.33: Interferogram between two bistatic radar images [109]. 


measurements. This means that the up and down-frequency converters did 

not introduce uncontrolled phase shift. 
It is important to note that the same values could not be obtained with- 
out the calibration path provided by the third antenna of the transponder. 
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Figure 3.34: Time series of the three measured displacements: monostatic, 
bistatic without calibration, and bistatic with calibration [109]. 


Indeed, the bistatic radar image was acceptable even without calibration of 
the transponder, but the mean value of the displacement of point A resulted 
in mean value of 2.840 mm with a standard deviation of 4.247 mm, as shown 
in Figure 3.34. These values prevented any possible use of the interferograms 
for detecting displacements. 

The reason for this phase instability is related to the coaxial cable that 
connects the radar head to the transponder for sharing the 10 MHz reference. 
When the radar head moves along the mechanical guide, the movement of 
the cable gives uncontrolled phase changes. 


Cross-pol long-cable Transponder 


An other possible solution to increase the isolation between the two antennas 
of the transponder is to enlarge the mutual distance. This solution is suitable 
also when the line of sight between radar and transponder is not present. 
In this case the two antennas are cross-polarised, as well as shown in Fig- 
ure 3.19 and 3.35. An RF cable separates the two antennas of the transpon- 
der. The transponder has two amplifiers A, and A» at the two ends of the 
cable. The gain of A; is 84 dB and the gain of Ao is 28 dB. The loss of the 
RF cable is - 30 dB. Therefore, the effective transponder gain is 82 dB. 
Figure 3.36 shows the aerial picture of the test site, where the radar 
equipment was tested. The radar was facing a seven-storey building at 140 
m distance. The TX antenna of the transponder was at the right-hand side 
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Figure 3.35: Working principle of monostatic/bistatic GBSAR using cross- 
polarised long-cable transponder. 


of the radar at 53 m distance. The RX antenna of the transponder was 
at 31 m distance from the TX antenna. The radar operated in the band 
9.915-10.075 GHz. The number of frequencies was 801. The length of the 
mechanical scan (orthogonal to the view direction) was 1.82 m. The number 
of steps along the scan was 180. The transmit power was 19 dBm. 


Figure 3.36: Aerial picture of test site. 


Figure 3.37(a) shows the obtained monostatic image. The edges of the 
building are well-evident. The bistatic image is shown in Figure 3.37(b). 
The arrow indicates the same point P in the aerial picture, the monostatic 
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image, and the bistatic image. 
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(a) Monostatic radar image of test site (b) Bistatic radar image of test site [113]. 
[113]. 


Figure 3.37: Monostatic/bistatic radar image of test site [113]. 


The TX antenna of the transponder and the monostatic RX antenna 
of the radar are cross-polarised, so their coupling should be low. Never- 
theless, the transponder (with 82 dB gain) is active even when the radar 
operates in monostatic modality. Therefore, a week signal can pass through 
the transponder producing some artefacts. To verify it, the bistatic signal 
has been focused with the monostatic algorithm. The obtained image is 
shown in Figure 3.38. The feature marked with A in the monostatic image 
(Figure 3.37(a)) can be the “ghost” image of the building. Its amplitude is 
low (40 dB lower than the highest signal), but it must not be mistaken for 


a true target. 
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Figure 3.38: Bistatic signal focused with monostatic algorithm [113]. 


To test the phase stability of the equipment, 12 images (monostatic and 
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bistatic) in 8 h were acquired and the interferograms between any image and 
its subsequent image has been calculated. As an example, Figure 3.39 shows 
one monostatic interferogram. The pixels relative to the building appear 
rather stable in phase. 
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Figure 3.39: Example of monostatic interferogram [113]. 


To quantify this phase stability, the “apparent” displacement of the point 
marked with P has been calculated as (2.22). 

As the point marked with P is presumably a stable point, any displace- 
ment fluctuation is due to the radar system, atmospheric fluctuation, and 
environmental temperature changes. Figure 3.40 shows the plots of these 
displacements for monostatic and bistatic configurations. These statistical 
fluctuations are in good agreement with other measurement campaigns with 
similar equipment [7], [112]. 
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Figure 3.40: Displacement of point P as detected by 11 monostatic and 
bistatic interferograms [113]. 
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3.5 Other applications 
Bistatic ArcSAR 


ArcSAR is a particular implementation of GBSAR, which has been receiving 
increasing interest in the last years. The great advantage of ArcSAR is its 
capability to synthesize images at 360° with constant resolution in azimuth 
[24]. 

Figure 3.41 shows a sketch of the radar prototype assembled for testing 
the working principle of bistatic ArcSAR [105]. The vector network ana- 
lyzer (VNA) HP8720D operated as transceiver providing a continuous wave 
stepped frequency signal (SFCW) in X-band with central frequency fe = 
10 GHz and bandwidth B =160 MHz. Two RF cables linked the VNA to 
the front-end fixed at the tip of the rotating arm. 

Before the receiving chain a switcher selected the signal from the mono- 
static/bistatic antenna. The bistatic antenna (omnidirectional) was con- 
nected to the front-end trough a 25 m long RF cable. Two 40 dB amplifiers 
compensate the cable loss. 


RADAR 
FRONT-END 


OMNIDIRECTIONAL 
ANTENNA 


` 
` 
SWITCHES’, 
CONTROL * 
i 


Figure 3.41: Bistatic ArcSAR prototype [105]. 


For testing the capability to obtain monostatic and bistatic radar images 
the ArcSAR and the transponder were installed in a garden of the University 
of Florence. In front and on the side of the radar there was a metallic fence 
as reported in the map in Figure 3.42. 

Figure 3.43 shows the obtained monostatic and bistatic images. 

Both images appear of good quality. As expected the resolution and the 
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Figure 3.42: Map of the radar installation [105]. 


signal-to-noise ratio of the monostaticimage are better. 
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(a) Monostatic image [105]. 
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(b) Bistatic image [105]. 


Figure 3.43: Monostatic/Bistatic ArcSAR image. 
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With the aim to proof the capability of the bistatic radar to detect two 
components of the displacement a special target was assembled [105]. The 
target was mounted over a linear micrometric positioner for optical bench 
with 0.1 mm nominal accuracy. 

The positioner was aligned along the y-axis. A reference couple of im- 
ages (monostatic and bistatic) was acquired. Aftermath, the positioner was 
moved of 10 mm along y-axis, and a second couple of images was acquired. 


The obtained displacements along y and x were +1.1 mm and -7.61 mm 
(the nominal values were 0.0 mm and -10 mm). The agreement is not perfect 
but the working principle has been demonstrated also in this case. 


Dynamic measurement for monitoring slander structures 


The bistatic approach can be used also for monitoring slander structures 
like telecommunications towers, wind turbine towers, bell-towers, minarets, 
obelisks and other monuments. Indeed this slander vertical structures need 
both testing before going into service and periodic health monitoring. Their 
measured responses when subjected to dynamic stimula (e.g. wind or vehic- 
ular traffic) is a key parameter for assessing their structural integrity [114]. 

By using a bistatic approach it is possible to detect, at the same time, 
two components of the dynamic displacement and it has been possible to 
identify the axis of symmetry of the structure. Any rotation of this axis 
can be associated to a possible structural change and can be considered a 
warning for deeper investigations [115]. 

Figure 3.44 sketches the working principle of the monostatic/bistatic 
radar. An interferometric radar positioned in A detects the component of 
the displacement along the direction target-A (for any point scatterer of the 
tower). A transponder is positioned in B. The radar in A transmits a signal 
that is scattered by the tower, bounces through the transponder and comes 
back to the radar. As the electromagnetic wave is transmitted from A point 
and received in B point, the radar detects the component of the displace- 
ment at the bistatic angle (i.e. the bisector of the angle A-target-B). As the 
radar is able to perform the two measurements nearly simultaneously (4.3 
ms of delay), the effective displacement vector in time (for any point target 
of the tower) is retrieved. 

The radar is a modified version of IBIS-FM (see Figure 3.45), produced 
by IDS GeoRadar, Pisa, Italy [116]. It operates an FM signal in Ky band- 
width, with maximum operative bandwidth of 400 Hz. The radar is provided 
with 2 transmitting (TX) and 2 receiving (RX) channels. The radar can op- 
erate in standard configuration, by using one TX and one RX antenna, or 
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Figure 3.44: Working principle of monostatic/bistatic radar technique ap- 
plied to the dynamic monitoring of a tower. 


in MIMO (Multiple Input Multiple Output) configuration by using all the 
combinations of the four antennas sequentially. The transponder consists of 
a RX-antenna, a large-bandwidth amplifier (Gain: 56 dB), a 3.0 m RF-cable 
(attenuation: -5 dB), and a TX-antenna. The two antennas can be pointed 
using two tripods. 

The radar system has been tested on a telecommunication tower in Flo- 
rence, Italy. It is steel tower of 35 m height. Its basement has 1.17 m 
diameter. Figure 3.46 shows the measurement setup: on the left of the 
transponder, on the right the radar, and in the background the tower. A 
light wind was the only mechanical stimulus. 

The radar operated in MIMO modality. The sampling frequency for 
each channel was 234.7 Hz. The starting frequency of radar was 17 GHz, 
the stopping frequency was 17.4 GHz with 1051 time-samples. The range 
resolution of the equipment was about 0.38 m. The set-up geometry is shown 
in Figure 3.47. The bistatic angle was 27.37°. The distance between radar 
and tower’s basement was 9.82 m. 

For obtaining the radar plot in time the Fast Fourier Transform (FFT) 
of the radar echo was calculated. In order to retrieve the displacement of 
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Figure 3.45: Modified version of IBIS-FM. 


Figure 3.46: Measurement configuration. 


the same target, the monostatic and bistatic plots have been projected on 
the tower axis by considering the electromagnetic phase history (see (2.18)). 

With reference to Figure 3.48, the path between a generic k-point on the 
vertical axis is dj/°"° = 2Rq,,) for the monostatic, and dP‘? = Rak) + 
Rear) + Leable + Rr, for the bistatic. 

Figure 3.49 shows the plots obtained by compensating the paths. They 
show a notable overlapping, but a perfect correspondence of all peaks is not 
possible, as each target could have a different response in the monostatic 
and bistatic configuration. The target at 19.55 m height has been selected 


70 


LAPO MICCINESI 


s-O” N Tower 


Transponder - R; 


RF-Cable 


2.91m 


Transponder - TX 


afr? = 2- Ray 


dpi = Rye + Roe + Leapte+ Rr 


{1 i] | | jj LILLLILII 
= HEATH HHHH 


. Rr 
Leable 


Rx antenna Tx antenna 


Figure 3.48: Measurement configuration. 


for retrieving the displacement. 

The duration of the whole measurement was about 15 minutes. Figure 
3.50 is an example of only 18 s. We can observe that the two displacements, 
monostatic and bistatic, are in phase. 

For retrieving the natural frequency, the Joint Time Frequency Analysis 
(JTFA) was applied to the displacement [117]. The Fourier transform of 
displacement has been evaluated in a sliding-windows of 76.83 s with an 
overlapping of 90%. 

The spectrum of natural frequency (calculated as average in time of the 
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Figure 3.49: Radar plot for monostatic (blue) and bistatic (red). 
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Figure 3.50: Comparison bewtween displacement measured by mostatic 
(blue) and bistatic (red) in a period of 18 s. 


absolute value of JTFA) is shown in Figure 3.51. The natural frequency of 
the tower was in 0.65 Hz both for monostatic and bistatic. 


The x and y components of displacement are shown in Figure 3.52. 


In order to find the axis of symmetry of the tower, the x-y components 
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Figure 3.51: Natural frequency of the tower obtained by a time-avarage of 
the amplitude of the JTFA (mostatic-blue, bistatic-red). 
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Figure 3.52: x (blue) and y (red) components of displacement of the tower. 


have been reported in a scatter plot by considering the full acquisition time. 
A liner fit has been applied on the scatter plot in order to estimate the 
direction of axis of symmetry (Figure 3.53(a)). The axis of symmetry can 
be estimated as (-17.59 + 0.05)°. 
A better estimation of axis of symmetry can be obtained by filtering the 
signal at the natural frequency using a bandpass Butterworth filter with 6 
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poles with central frequency 0.65 Hz and band 0.04 Hz. The axis of symmetry 
of the filtered displacement is shown in Figure 3.53(b). The points stand now 
precisely along a straight line. So a better evaluation of the axis results to 
be (- 16.218 + 0.012)°. This direction corresponds to 102.8° North. 
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Figure 3.53: Axis of symmetry of the tower considering whole bandwidth 
(a) and filtered (b). 


3.6 Conclusions 


In this chapter a GBSAR operating in monostatic and bistatic modality has 
been presented. A transponder was made use in order to obtain a secondary 
image, and one other component of displacement, from different point of 
view. 

The radar has been tested with a basic transponder (two antennas and 
an amplifier) in a controlled environment and it was able to detect two 
components of displacement. 

In order to improve the SNR of the bistatic measurements two enhanced 
transponders, with an hight gain, were presented. 

The first transponder made use of a frequency shifter in order to insulate 
the two antennas and achieved a gain of 91 dB. The artefacts did not affect 
the measurement with this transponder, but its implementation is rather 
difficult in practical cases. 

The other is a long cable transponder with a gain of 82 dB. This transpon- 
der appear to be more practical for most of the uses. 

Other application of bistatic technique was presented. A notable appli- 
cation is monitoring slander structures. In this application it was possible 
to measure the natural axis of a telecommunication tower. 
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Compressive Sensing MIMO Radar 


The working principle of GBSAR is based on dense spatial sampling (steps 
smaller than a quarter of wavelength) along the synthetic aperture (usually 
a linear mechanical guide) as explained in Chapter 1. 

Compressive Sensing (CS) is a recent sampling paradigm [118], [119] 
which asserts that it is possible to recover certain signals from far fewer 
samples or measurements than traditional methods use. Its basic idea relies 
on the “sparsity” of the signals of interest (the radar signals typically have 
this property [120]), and the incoherence of the sensing modality. The latter 
property is obtained through random sampling. 

Generally speaking, CS can be applied in frequency domain (for reducing 
the number of transmitted frequencies) or in spatial domain (for reducing 
the number of spatial steps along the synthetic aperture). 

In 2010, Huang et al. [121] applied CS to step-frequency through-wall 
radar that scanned along a horizontal mechanical guide. In 2011, Karlina and 
Sato [122] proposed to apply CS to GBSAR and they performed simulations, 
as well as Zonno in 2014 [123]. Yigit [124] extended CS to GB-SAR/ISAR ac- 
quisition modality and performed measurements in controlled environment. 
The interferometric properties of CS-GBSAR have been demonstrated by 
Giordano el al. [125] in 2015. 
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At the state of the art, CS techniques appear very promising for GBSAR 
applications, but all the cited works are relative to radar acquisitions in 
controlled environments (anechoic chamber or short-range experiments with 
controlled targets). Therefore, the crucial question is if these techniques 
are effectively applicable with data acquired in real scenarios, like GBSAR 
images of large buildings, landslides, open-pit mines, and glaciers. 

Furthermore, a few specific questions about the CS methods have to be 
addressed. The first one is how to choose the orthogonal set of functions 
(the basis) that are used for recovering the signal. The possible bases are 
countless. The possible recovery algorithms are many as well. Therefore, a 
general selection criterion is necessary both for basis selection and recovery 
method. Another question is: what is an acceptable compression for real 
data. The claims declared in works that use simulated data or acquisitions 
performed in very cooperative scenarios are probably too optimistic. 

Finally, in previous papers [123], [125] it is possible to note that the y 
matrix (see next paragraph for his definition) is built by under-sampling the 
complete set of data acquired along a single scan. It means that each line 
could use a different set of positions along the scan. In this article, the y 
matrix is built using the same reduced set of data for all the lines. This could 
reduce the incoherence of sampling (and so the quality of reconstruction), 
but it is a solution that could be easily implemented with a physical array or 
with a MIMO (multiple input multiple output), such as in [126], [25], [26]. 

In this Chapter the theory of Compressive Sensing will be introduced. 
After, a systematic study about basis, reconstruction algorithm and the 
acceptable compression value will be presented [127]. Finally, a Compressive 
Sensing Interferometric MIMO Radar is introduced [30] with some of its 
application for bridges monitoring [128]. 


4.1 Compressive sensing 


Compressive sensing theory states a sparse signal can be recovered by fewer 
samples than what is required by the Nyquist theorem. With reference to 
Figure 4.1, a standard GBSAR transmits and receives at N positions along a 
linear mechanical guide. The Nyquist theorem requires that the spatial step 
has to be smaller than a quarter of wavelength (A/4) for omnidirectional 
antennas (this constraint is a bit more relaxed for directional antennas, but 
it is not essential in the discussion that follows). Let E be the vector of the 
N samples acquired according the Nyquist theorem. A subset of M random 
positions is selected between the N positions. These M positions are the 
effective points where the radar-head performs the measurements that will 
be used for recovering the complete set of N samples. 
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Figure 4.1: Sampling along the linear mechanical guide of a GBSAR. 


The next step is to define M variables (ym ) as random linear combinations 
of the N samples using only the M positions as sketched in Figure 4.2. 
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Figure 4.2: Definition of the M variables y,,. The number different from 
zero in the matrix are selected randomly between 0 and 1. 


The vector y of the M variables can be written as 
y= gE (4.1) 


with y matrix M x N as defined in Figure 4.2. It should be noted that 
the vector y is a measurable quantity, as it can be calculated from the M 
measurements. 

A radar signal is often sparse in a suitable set of functions. As an ex- 
ample, an ideal scenario with a single point-target in far field appears as a 
single peak in the domain of spatial frequencies. Inspired by this idea, we 
can write the vector E as 


E = qb (4.2) 
77 


ADVANCED GROUND-BASED REAL AND SYNTHETIC APERTURE RADAR 


where w can be 


1 0 0 


and where b is the vector of the coefficients of spatial Fast Fourier Transform 
(FFT) of the signal. If the signal is sparse, many coefficient of basis b, are 
almost null. 

The FFT or the DCT (discrete cosine transform) should be able to ef- 
fectively represent ideal point targets, but they could be a too rough ap- 
proximation for more realistic distributed targets. Wavelets transforms give 
better results with distributed targets [129]. The possible complete sets of 
wavelets are countless, so this study has been limited to the most popular 
sets [129], [130], [131], [132] specifically: Haar wavelets, Daubechies wavelets, 
biorthogonal wavelet, coiflets, symlets, LeGall wavelets and discrete Mayer 
wavelets. 

By substituting (4.2), in (4.1), the vector y can be expressed as 


y = vb (4.4) 


with J = pw. This matrix is the so-called “dictionary”. 

Equation (4.4) is a linear system with many more variables than equa- 
tions. It is an ill-posed system but it can be solved with a suitable recovery 
method. The most popular are: 11-minimization [133], 12-minimization [134], 
and orthogonal matching pursuit [135]. Using one of these methods, the vec- 
tor b is recovered, so the E vector of N samples is obtained by (4.2). 

The procedure described above has to be repeated for each frequency for 
obtaining the matrix Ez, with k-index ranging from 1 to Ny and i-index 
ranging from 1 to N. 

The next step is to focus the matrix Ex, using the back-propagation 
algorithm described in Chapter 2. 


4.2 Compressive sensing for GBSAR 


In this section a systematic investigation of the possible dictionaries and 
recovery algorithms is presented. 

In order to estimate the quality of the recovered images the peak signal- 
noise ratio (PSNR) was calculated. The PSNR is defined as 


max(|I;,;|2) 


PSNR = 10logio a 
a oS dag ag? 
t=1j=1 


(4.5) 
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where R; j is the focused reference GBSAR image obtained from the data 
sampled according to Nyquist theorem without application of CS techniques; 
I; j is the image obtained using the CS techniques. Both images are complex 
matrix m x n. As rule of thumb, images with PSNR lower than 20- 25 dB 
appear of low quality, with features hardly recognizable [136]. 


Simulations 


First of all the minimum compression has to be founded. For this reason, 
the CS algorithms have been preliminary tested with the simulation of a 
single target at 150 m distance in front of radar. The used parameters were: 
initial frequency fı = 9.915 GHz, final frequency fa = 10.075 GHz, number 
of frequencies N; = 801, length of the mechanical scan L = 1 m, number of 
samples (according to Nyquist theorem) N = 100. 

The basis used was the Haar wavelets. The recovery method was 12- 
minimization (this choice of Haar wavelets will be motivated in the experi- 
mental section, 12-minimization has been selected because it works even with 
a number of variable smaller than the number of equations, nevertheless the 
choice of both basis and recovery method is not essential for the discussion 
that follows). The PSNR has been calculated for an image 200 x 200 pixels 
in polar coordinates. The range limits were 100 m and 200 m. The azimuth 
limits were -30 deg and +30 deg. 

Figure 4.3 shows the plots of PSNR versus M/N. Each PSNR was cal- 
culated as the average of 25 random iterations and the standard deviation 
of each value is reported in the plot (blue line). It is interesting to note 
that for M/N > 0.3, the PSNR increases linearly of 2.4 dB for each 10% 
increase of M/N. This is a notable result that is also true for simulations 
of more complex scenarios (several targets), using different measurement 
parameters, different dictionaries, and different recovery methods, and in 
experimental data too. Note that PSNR does not tend to infinity for M = 
N, but it stops at about 37 dB. The reason of it is related to the limited size 
of the y matrix. 

As it is known, the CS reconstruction relies on the incoherence of both 
random sensing and the random coefficients of the matrix y. When M = N, 
the sensing is no longer incoherent, but the incoherence of the y coefficients 
persists and limits the quality of reconstruction. The y coefficients are uni- 
formly distributed random numbers between 0 and 1. When M = N, the 
number of elements is N x N = 100 x 100 = 10000, the sum of 10000 uni- 
formly distributed random numbers between 0 and 1 gives exactly 37.0 dB, 
that is just the maximum PSNR obtained in Figure 4.3. 

The noise can affect both PSNR and displacement retrieved by inter- 
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Figure 4.3: PSNR in function of M/N (simulated data) [127]. 


ferometry. In order to verify/evaluate it, we have repeated the simulations 
reported above with controlled Gaussian noise. Figure 4.4 shows the PSNR 
versus M/N calculated with different SNR values (the SNR reported in the 
plots is relative to one single measurement specified by one single frequency 
and one single position along the mechanical guide). Each PSNR value of 
the plots in Figure 4.4 is calculated as the average of 12 random iterations. 
As expected, by reducing the SNR the PSNR decreases. When the SNR of 
a single measurement is > 0, the PSNR is practically equal to the PSNR 
without noise calculated in the previously. 


a SNR = 0dB 
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Figure 4.4: PSNR in function of M/N and SNR (simulated data) [127]. 


The noise can give an error in the displacement retrieved by interfer- 
ometry. In order to verify/evaluate it, we have estimated the phase error 
(A¢@) and we have converted it in displacement error (AR) using the well- 
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known relationship (2.22). 


Figure 4.5 shows the estimated error in function of SNR. Each PSNR 
value is calculated as the average of 40 random iterations. For N4/N = 1 
the error plot in function of SNR overlaps perfectly the plot of the error 
calculated without CS that, in log scale, is a straight line with slope -1. As 
expected, for M/N = 0.5 the error increases of 3dB and for M/N = 0.1 the 
error increases of 10 dB. 
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Figure 4.5: Displacement error in function of SNR for different M/N (simu- 
lated data). The dark line is the displacement error calculated without CS 
technique [127]. 


When M/N increases, the number of samples increases too. So the error 
should decrease. This intuitive idea has been verified with a simulation. 
Figure 4.6(a) shows how the error decreases increasing M/N, for two different 
values of the SNR. Both the plots are approximately straight lines with 
slope -1. 


Generally speaking, the interferometric error could not depend only on 
Gaussian noise (i.e. on thermal noise). In that case, calculating the error in 
function of M/N does not make much sense, as it could be hard to obtain 
a straight line. On the contrary, the differential error (i.e., the error with 
respect to the image obtained using the CS with 100% of the samples) is 
expected to have a linear behaviour (in log scale). Figure 4.6(b) shows the 
differential displacement in function of M/N for three different SNR. It is 
evident that the plots can be approximated as straight lines with a slope 
of -2. 
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(a) Displacement error in function of M/N for dif- 
ferent SNR (simulated data) [127]. 
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(b) Differential displacement error in function of 
M/N for different SNR (simulated data) [127]. 


Figure 4.6: Displacement error and differential displacement error in function 
of M/N for different SNR (simulated data) [127]. 


Experimental application 
Sesto Fiorentino Campus 


The CS techniques have been experimentally tested in a controlled scenario 
at the University of Florence (named Sesto Fiorentino Campus). A single 
metallic target, provided with a micrometric linear translator, has been po- 
sitioned in front of the radar equipment at 150 m distance in a large flat 
garden, as shown in the pictures in Figure 4.7. The radar operates in mono- 
static modality. 

The measurement parameters were the same used in the simulations. 
The acquired experimental data have been used for testing the most popular 
bases and recovery methods. Specifically we used discrete cosine transform 
(det), fast Fourier transform (fft), the Haar wavelts (haar), the Daubechies 
wavelets from order 2 to order 10 (db2-db10), seven biorthogonal wavelets 
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Figure 4.7: Experimental test site Sesto Fiorentino Campus [127]. 


(bior1.3, bior2.2, bior2.6, bior3.1, bior3.5, bior3.9, bior5.5), the “coiflets” 
from order 1 to 5 (coifl-coif5), the “simlets” from order 2 to 8 (sim2-sim8), 
LeGall wavelets (legal5.3), the discrete Meyer wavelets (dmey). As recovery 
methods, we used ]1-minimization (L1), 12-minimization (L2), and orthogo- 
nal matching pursuit (OMP). 

Table 4.1 reports the obtained PSNR for M/N = 0.5 and M/N = 0.33. 
The PSNR has been calculated for an image 200 x 200 pixels in polar coor- 
dinates. The range limits were 100 m and 200 m. The azimuth limits were 
-30° and +30°. 

In the table, the highest (in bold) and the lowest (in italic bold) values 
were highlighted. Using 50% of data, PSNR varies between 26.01 dB (Haar 
wavelets, OMP) and 21.89 dB (FFT, OMP). Using 30% of data, PSNR varies 
between 20.03 dB (Haar wavelets, OMP) and 15.87 dB (Bior3.1 wavelets, 
OMP). 

In simulated data, it is possible to note that PSNR increases linearly of 
2.4 dB for each 10% increase of M/N. This relationship is perfectly confirmed 
also for these experimental data as shown in Figure 4.8. 

For evaluating the interferometric error 100 measurements in 96 h were 
acquired. By considering one pixel associated to a stable point of one build- 
ing at about 200 m in the background, the standard deviation of the dis- 
placement retrieved by interferometry has been evaluated. Figure 4.9 shows 
the obtained values varying M/N. The error is rather constant. It means 
that thermal (Gaussian) noise is not the main source of error. The effective 
contribution of thermal noise is more evident in the plot of the differential 
error (Figure 4.10) that confirms the error (in log scale) linearly decreases 
with slope -2, increasing the M/N (in log scale), as predicted by simulation. 
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50% 30% 
L1 (dB) L2(dB) OMP (dB) L1(dB) L2(dB) OMP (dB) 

Det 22.92 22.89 22.29 18.77 19.02 19.62 
Fft 22.61 22.73 21.89 17 17.16 17.45 
Haar 24.83 24.83 26.01 18.3 18.3 20.03 
Db2 24.32 24.29 24.32 18.14 18.13 18.13 
Db3 24.32 24.33 24.36 17.68 17.66 18.01 
Db4 23.97 23.97 23.67 18.22 18.23 17.98 
Db5 23.71 23.73 23.1 17.49 17.48 17.35 
Db6 23.87 23.9 22.74 17.76 17.77 17.98 
Db7 23.24 23.25 22.6 17.41 17.45 16.11 
Db8 23.24 23.23 22.48 17.71 17.7 17.82 
Db9 23.26 23.27 21.46 17.71 17.72 17.82 
Db10 23.04 23.07 23.48 17.74 17.74 18.52 
Bior1.3 24.65 24.64 25 18.16 18.16 19.51 
Bior2.2 24.45 24.45 23.53 17 16.99 16.96 
Bior2.6 24.22 2421 23.37 17.65 17.65 17.37 
Bior3.1 23.51 23.51 22.66 17.39 17.39 15.87 
Bior3.5 24.29 24.28 23.5 17.99 17.99 17.21 
Bior3.9 24 24.04 23.67 17.6 17.6 17.67 
Bior5.5 23.97 23.98 23.62 18.04 18.07 18.23 
Coifl 24.04 24.09 23.52 17.71 17.72 17.25 
Coif2 24.2 24.21 24.04 16.84 16.84 17.26 
Coif3 23.72 23.72 23.27 17.39 17.39 17.36 
Coif4 23.47 2347 22.74 17.49 17.51 17.52 
Coif5 23.28 23.18 22.79 17.05 17.04 16.9 
Sym2 24.32 24.29 24.32 18.14 18.13 18.13 
Sym3 24.32 24.33 24.36 17.68 17.66 18.01 
Sym4 24.42 24.44 24.08 17.29 17.28 17.04 
Sym5 23.77 23.77 23.62 17.63 17.64 17.54 
Sym6 24.24 2415 23.78 17.64 17.64 16.94 
Sym7 24.26 2427 24.18 18.02 18.03 17.68 
Sym8 24.1 24.11 23.53 17.94 17.94 17.04 
Legall5.3 2417 24.17 24.49 18.31 18.3 18.51 
Dmey 23.6 23.6 23.56 17.93 17.94 17.96 


Table 4.1: PSNR calculated with different bases and recovery methods for 
experimental data acquired at the test site Sesto Fiorentino Campus. 
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Figure 4.8: PSNR in function of M/N (data acquired at the experimental 
test site Sesto Fiorentino Campus) [127]. 
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Figure 4.9: Displacement error in function of M/N for a stable point in the 
test site Sesto Fiorentino Campus [127]. 


The metallic target was provided with a micrometric positioner able to 
translate it with a nominal accuracy of 0.1 mm. So, we performed 20 radar 
acquisitions by moving the target of 1 mm forward the radar after each 
acquisition. Figure 4.11 shows the cumulative displacement detected by 
interferometry using 100% of data without CS, using 100% of data with CS, 
using 50% of data with CS, and using 30% of data with CS. The four plots 
are practically overlapping. Figure 4.12 shows the plots of the differences. 
In the worst case (CS with 30% of data) the error is lower than 25 um. 
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Figure 4.10: Differential displacement error in function of M/N for a stable 
point in the test site Sesto Fiorentino Campus [127]. 
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Figure 4.11: Cumulative measured displacement of the controlled metallic 
target at the test site Sesto Fiorentino Campus [127]. 


A Seven-Storey Building 


In order to test the application of CS techniques in an operative scenario 
representative of an urban environment, we performed radar measurements 
using as target the seven-storey building shown in Figure 4.13. 

The measurement parameters were: initial frequency fı = 9.915 GHz, 
final frequency f2 = 10.075 GHz, number of frequencies N; = 801, length 
of the mechanical scan L = 1.8 m, number of samples (according to Nyquist 
theorem) N = 180. The obtained radar image (without application of CS 
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Figure 4.12: Differences between displacements detect through 100% CS, 
50% CS, 30% and displacements detected using 100% of data without appli- 
cation of CS techniques (test site Sesto Fiorentino Campus) [127]. 


Figure 4.13: Aerial picture of the seven-storey building used as urban test 
site [127]. 


techniques) is shown in Figure 4.14 


As well as for simulated data, and for the experimental data acquired at 
the “Sesto Fiorentino Campus” test site, the PSNR increases linearly of 2.4 
dB for each 10% increase of M/N and this relationship has been perfectly 
confirmed also for these experimental data (for sake of brevity the graph is 
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Figure 4.14: Radar image (without application of CS techniques) of the 
building shown in Figure 4.13 [127]. 


not reported). 

The PSNR for M/N = 0.5 and M/N = 0.3 has been calculated (Table 
2 in [127]). Using 50% of data, PSNR varies from 22.78 dB (db2 wavelets, 
OMP) to 8.19 dB (dct, OMP). Using 30% of data, the PSNR ranges between 
19.12 dB (db5 wavelets or sym2 wavelets, OMP) and 1.80 dB (fft, OMP). 


An Open-Pit Copper Mine 


The most popular use of the GBSAR systems is the monitoring of open- 
pit mines. Therefore, the CS techniques has been tested to a representative 
case study of a copper mine in South America. The measurement parameters 
were: initial frequency fi = 17.1 GHz, final frequency fa = 17.3 GHz, number 
of frequencies N; = 5333, length of the mechanical scan L = 1.275 m, number 
of samples (according to Nyquist theorem) N = 256. The obtained radar 
image (without application of CS techniques) is shown in Figure 4.15. 
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Figure 4.15: Radar image (without application of CS techniques) of a copper 
mine in South America [127]. 
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As well as for previous data, the PSNR increases linearly of 2.4 dB for 
each 10% increase of M/N and this relationship has been perfectly con- 
firmed also for these experimental data (for sake of brevity the graph is not 
reported). 

Also in this case the PSNR for M/N = 0.5 and M/N = 0.3 has been 
evaluated (Table 3 in [127]). Using 50% of data, PSNR varies from 29.95 
dB (Haar wavelets, OMP) to 17.25 dB (FFT, L2). Using 30%, the PSNR 
ranges between 29.44 dB (Haar wavelets, OMP) and 9.33 dB (FFT, OMP). 

Figure 4.16 shows an interferogram between two images (focused with- 
out application of CS techniques) taken at different time. It shows a clear 
movement on the left side of the slope. In order to test the capability of 
the CS techniques to reconstruct interferograms, it is applied to 50% of data 
(Figure 4.17(a)) and to 33% of data (Figure 4.17(b)) using Haar wavelets 
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Figure 4.16: Interferometric image (without application of CS techniques) 


of a copper mine in South America [127]. 
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(a) Interferometric radar image obtained (b) Interferometric radar image obtained 
with 50% of data and CS [127]. with 33% of data and CS [127]. 


Figure 4.17: Interferometric radar image obtained with 50% and 33% of data 
and CS [127]. 
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In order to evaluate the interferometric error, 41 measurements acquired 
in 4 h and 20 min have been processed. By considering one pixel associated 
to one stable point (without visible fringes), the standard deviation of the 
displacement retrieved by interferometry has been evaluated. The point 
is labelled with letter A in Figure 4.15. Figure 4.18(a) shows the error 
(calculated as standard deviation) of the displacement at the stable point A 
varying M/N. The error is rather constant. It means that thermal (Gaussian) 
noise is not the main source of error. The effective contribution of thermal 
noise is more evident in the plot of the differential error (Figure 4.18(b)). 
It confirms that the error linearly decreases (in log scale) with slope -2, 
increasing the M/N (in log scale), as predicted by simulation. 
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(a) Displacement error in function of (b) Differential displacement error in 
M/N for a stable point in the open pit function of M/N for a stable point in the 
mine [127]. open pit mine [127]. 


Figure 4.18: Displacement and differential displacement error in function of 
M/N for a stable point in the open pit mine [127]. 


In order to verify the capability of CS to detect a cumulative displace- 
ment, we have taken into account the point B in Figure 4.15 where several 
fringes are evident in the interferograms (see Figure 4.16). Figure 4.19 shows 
the plot obtained without CS and the plot obtained with CS applied to 100% 
of data are perfectly overlapping. The CS plots with 50% and 30% of data 
tend to underestimate the displacement. 


A Glacier 


9 


As test-site representative of a natural scenario we considered the “Belvedere’ 
glacier in the Italian Alps. The data was acquired during a measurement 
campaign in 2006, whose details are reported in [137], [138], [139]. The 
measurement parameters were: fı = 5.97 GHz, f2 = 5.99 GHz, N; = 801, 
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Figure 4.19: Cumulative measured displacement of the point B in the open 
pit mine [127]. 


L = 1.71 m, N = 141. The obtained radar image (without application of CS 
techniques) is shown in Figure 4.20. 
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Figure 4.20: Radar image (without application of CS techniques) of a glacier 
in the Italian Alps [127]. 


As well as for previous data, we verified that PSNR increases linearly of 
2.4 dB for each 10% increase of M/N and this relationship has been perfectly 
confirmed also for these experimental data. 

The PSNR has been calculated for an image 200 x 200 pixels in polar 
coordinates. The range limits were 350 m and 1000 m. The azimuth limits 
were -25° and +25°. Using 50% of data, PSNR varies from 21.16 dB (Haar 
wavelets, OMP) to -3.50 dB (DCT, L2). Using 30%, the PSNR ranges 
between 16.52 dB (Haar wavelets, OMP) and -5.09 dB (FFT, L1). 

In order to verify the capability of CS to detect a cumulative displace- 
ment, we have taken into account a pixel inside the sliding body of the 
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Campus Building Mine Glacier 


fi [GHz] 9.915 9.915 171 5.97 
fo [GHz] 10.075 10.075 173 5.99 
N; 801 801 5333 801 
L [m] 1 18 1295. 171 
L [m] 1 18 12755 171 
N 100 180 256 141 


PSNR max 50% [dB] 26.01 22.78 29.95 21.16 


Table 4.2: Parameters of each measurement campaign and the highest value 
of PSNR obtained using 50% of data. 


glacier. Figure 4.21 shows the plots obtained without CS and with CS 
applied to 100%, 50%, and 30% of data. The plots are almost perfectly 
overlapping. 
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Figure 4.21: Radar image (without application of CS techniques) of a glacier 
in the Italian Alps [127]. 


Table 4.2 resumes the parameters of each measurement campaign and 
the highest value of PSNR obtained using 50% of data. 

The CS techniques are able to reconstruct GBSAR images of good quality 
(PSNR > 20-30 dB) with almost 30 — 40% of data. 

Nevertheless, the quality of reconstruction depends on the selected basis 
and recovery method. In a test site with just one single strong target, the 
best image has been obtained with FFT and OMP. In an urban scenario the 
best choice has been Daubechies wavelets and OMP. For an open-pit mine 
the highest quality image has been reconstructed with Haar wavelets and 
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OMP. For a natural slope (a glacier) the optimal choice has been again Haar 
wavelets and OMP. In all cases, Haar wavelets with OMP gave good results 
(even if not always the best). 

Therefore, the Haar wavelets with OMP is a reasonable choice for most 
practical cases. Furthermore, we found that CS reconstruction does not 
worsen significantly the interferometric error. The plots of the cumulative 
displacement recovered with 100%, 50%, and 30% of data overlaps almost 
completely the plot obtained without application of CS techniques. 

All these findings constitute the basic scientific starting point for design- 
ing optimal sparse array or multiple-input multiple-output (MIMO) radar 
able to exploit the CS techniques. 


4.3 Compressive Sensing MIMO Interferometric Radar 


The MIMO radar proposed in this Chapter exploits the random distribution 
of 4 transmitting (TX) antennas and 4 receiving (RX) antennas as shown in 
Figure 4.22. 
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Figure 4.22: Sampling along the linear mechanical guide of a GBSAR. 


When the TX antenna in the i-th position and the RX antenna the j- 
th position are both switched on (and all the others switched off), it is 
equivalent (in far field) to transmit and receive with one single “virtual” 
antenna along the median axis in the position (i+j)/2. Therefore for each 
combination of the TX and RX antennas a specific pattern along the median 
axis is defined. This pattern can be seen as a random sampling of the 
electromagnetic field backscattered by the targets in the field of view of the 
radar. 
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Selection of optimal pattern 


Generally speaking, Nrx antennas positioned along a guide with N/2 pos- 
sible positions and Nex antennas positioned along another guide with N/2 
positions, give theoretically possible configurations (N/2)NTx . (N/2)Nrx 
(patterns). 

Obviously many of these are repetitive and not practically feasible as it is 
not possible to put two antennas too close or in the same position. In order 
to select between them the patterns more suitable for the CS procedure we 
have established the following criteria: 


1. the antennas must keep their order. It means that each antenna is 
numbered and antenna number n has to be always on the left of an- 
tenna n+1; 


2. as the antennas have a physical size, between two nearby antennas 
there must be at least Ap empty positions; 


3. the number of “virtual” antennas along the median axis must be max- 
imum, i.e. Nrx x NRx; 


4. equal patterns must be eliminated. For doing it any single pattern 
along the median axis is identified with a binary number of N digit. The 
true-value corresponds to one position occupied by one virtual antenna. 
The digit numbers are ordered and the double numbers discharged; 


5. in order to exploit the whole synthetic aperture, the first and the last 
positions along the median axis must be occupied. 


As an example we consider Nrx = Nrx = 4, N = 40, Ap = 4. The 
compression factor results in 16/40= 0.40. The patterns that respect the 
six conditions declared above are 132. Figure 4.23 shows these 132 patterns 
as they appear along the median axis of the MIMO. The “virtual” antennas 
occupy the yellow positions while the blue positions are empty. 

All these patterns could work with CS techniques. Nevertheless, some 
simulations has been performed in order to select the best possible. The 
radar response of one single target at 50 m distance has been simulated. 

The used parameters were: initial frequency fı = 9.84 GHz, final fre- 
quency fo = 10.16 GHz, number of frequencies N; = 801. The obtained 
image (using the CS with Haar wavelets and OMP recovery method) has 
been compared with the image obtained sampling as required by Nyquist. 
The extremes of the image were from 25 m to 75 m in range and from -25° 
to 25° in azimuth. The comparison has been performed estimating the Peak 
Signal Noise Ratio (PSNR) calculated as (4.5). 
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Figure 4.23: Possible patterns of “virtual antennas” along the median axis 
of MIMO [140]. 


The black plot in Figure 4.24 shows the PSNR for the 132 patterns as 
identified in Figure 4.23. All patterns give PSNR > 24 dB and the maximum 
value is about 29.5 dB. 


The performances of CS reconstruction depends on the specific image, 
generally they are better when the target is perfectly in front of the radar. In 
order to consider a more realistic scenario, the simulation has been repeated 
by shifting the target from x = -6 m to x = 6 m at 2 m step. The coloured 
plots in Figure 4.24 shows the PSNR for these targets. 


For selecting the best pattern a relative threshold has been selected. The 
threshold was lower than 1, but close to 1. The patterns with PSNR higher 
than its highest value, multiplied the relative threshold, has been identified 
for each plot. After that, only the pattern that are over the threshold for 
all the plots has been selected. Finally the threshold was raised until obtain 
a single pattern. With reference to Figure 4.23 the selected pattern is the 
number 15. 
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Figure 4.24: Simulated PSNR for different patterns and different target 
positions [140]. 


Radar prototype 


Figure 4.25(a) shows the block scheme of the radar prototype. The vector 
network analyser (HP8720D) operates as transceiver. 

The switching system consists of eight single-pole double-through (SPDT) 
mod. MSP2T-18-12+ and ten high phase stability (+0.5°) microwave cables 
(SUCOFLEX 126). Another couple of SPDT provides a calibration path 
(with a -40 dB power attenuator). 


A relay board controls the switching system. The calibration path and all 
the paths between VNA and each antennas are of the same electromagnetic 
length in order to avoid any further calibration procedure. Figure 4.25(b) 
shows a photograph of the radar head. 

The MIMO radar has been preliminary tested in a controlled experimen- 
tal test site. In an open garden, a single corner reflector (CR) of 0.4 m side 
length was positioned at 13.4 m in front of the radar. 

The measurement parameters were: fı = 9.84 GHz, fo = 10.16 GHz, 
and Ny = 801. Figure 4.26 shows the obtained radar image. 

The red line in Figure 4.27(a) is the Point Spread Function (PSF) in 
azimuth at the distance of the CR. The blue line is the PSF obtained by 
simulating one point scatter in front of the radar. The agreement between 
the two plots is very good. 

As the CS performance depends critically on the view angle, we repeated 
the measurement moving the CR 2.80 m on the right. Figure 4.27(b) shows 
the measured and simulated power plots. Note that the angle £ is between 
the x-axis and the view direction, so 90° is in front of the radar and y < 90° 
is on the right side. 
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(a) Block scheme of the radar prototype. 


(b) Photograph of the radar head. 


Figure 4.25: Block scheme and photograph of the radar head. 
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Figure 4.26: Radar image of a CR in front of the radar [30]. 


It is interesting to note that the simulated and measured plots have a 
small misalignment that cannot be easily corrected. Indeed, it is an effect 
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(a) Target in the front of radar. (b) Target at 2.80 m on the right. 


Figure 4.27: Experimental (red line) and simulated (blue line) PSF in az- 
imuth of a CR in front of the radar [30]. 


of the non-linearity of the CS recovery. 

In order to verify this statement, a simulation has been performed: a 
single target at 13.4 m distance with yo azimuth angle varying from 80° to 
100° at step of 0.5°. The error dy in the azimuth recovered using the CS 
has the non-linear behaviour shown in Figure 4.28. 


80 82 84 86 88 90 92 94 96 98 100 
Po (deg) 


Figure 4.28: Error in the azimuth angle estimation by varying the azimuth 
angle [30]. 


A critical issue of any radar is the amplitude of the sidelobes. For a CS 
radar, this value could be rather erratic. 

With the aim of evaluate it, the response of a single target at 13.4 m 
distance at do azimuth was simulated. The plot in Figure 4.29 shows the 
peak to sidelobe ratio (P/SL) in function of go. 

Only the sidelobes inside the antenna lobe (+20°) have been evaluated. 
It results that if the target is inside a view cone +15° the P/SL is always 
larger than 10 dB. 
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Figure 4.29: P/SL for a single target at 13.4 m distance at the azimuth angle 
between 0° and 18° at 1° step [30]. 


In order to test the capability of MIMO to detect displacement by inter- 
ferometry, the CR was located on a micropositioner with 0.1 mm nominal 
accuracy. The CR was moved forward the radar at step of 5 mm. 

Figure 4.30 shows the obtained results. The agreement with nominal 
values is better than 0.2 mm. 
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Figure 4.30: CR displacements detected by radar [30]. 


In-field tests 


The radar prototype has been in-field operated for monitoring a pedestrian 
bridge in Poggibonsi, Italy and Vespucci bridge in Florence, Italy. In the 
second case the results of radar prototype has been compared with a real- 
aperture radar in Ky bandwidth. 
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Pedestrian bridge 


The pedestrian bridge is located in Poggibonsi, Italy. The main span is 47.5 
m long (Figure 4.31) and it is built with similar block of 2.5 m long and 
2.8 m wide. In 2007, this bridge was tested with an interferometric radar. 
The experimental results were published in [111] and [141]. 
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Figure 4.31: Pedestrian bridge of Poggibonsi, Italy. 


The MIMO radar was installed close to one of the two pillars of bridge, as 
shown in Figure 4.31 and Figure 4.32. The measurement parameters were: 
initial frequency fı = 9.84 GHz, final frequency f2 = 10.16 GHz, number of 
frequencies N; = 801. Therefore, the range resolution was 0.47 m and the 
unambiguous range 375 m. The antenna aperture in the horizontal plan was 
about +20°. The transmitted power was 12 dBm. The radar completed a 
single acquisition in 31.4 s. 


Figure 4.32: Photograph of the radar installation. 
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Figure 4.33 shows the obtained power image of the lower deck. The first 
nine transversal beams are well visible. The shape of the deck is recognizable. 
Unfortunately, due to the low transmitted power (12 dBm), the transversal 
beams are hardly visible in the portion of the deck farther than 24 m. 
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Figure 4.33: Radar image of the lower deck of the bridge [30]. 


The bridge was loaded with a small car (900 kg). It slowly went to 
the median point and came back. In order to evidence the effects of an 
asymmetric load the car was driven as on the left as it was possible (the 
bridge is wide 2.8 m, while the car was 2.6 m wide) as shown in Figure 4.34. 


Figure 4.34: Test of the bridge using a small car as load. 


Figure 4.35 shows the detected displacements (projected in vertical di- 
rection) of the points A (on the left side of the deck) and B (on the right side 
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of the deck) during the loading and downloading operations. In the x-axis, 
there is the number of acquisition (the radar acquired an image each 31.4 s) 
Both the points are at 16.18 m range, very close to the center of the span. 

The car was slowly driven until to the center of span and it stayed there 
for about 12 min. As expected the point on the left had a larger displacement 
(about 2.0 mm). The difference of displacement between left and right sides 
is about 1 mm. 

We note at the beginning of measurement session a slight uplift of the 
deck. In effect, as the pillars were about 10 m inside the riverbed (see Figure 
4.31), when the car entered the bridge from one side the center of the span 
rose. Furthermore, the car passed the center and came back, so the plot in 
time shows two peaks (at 21st and 33rd measurements). 


displacement along z (mm) 
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Figure 4.35: Plots of the displacements of points A and B [30]. 


Vespucci bridge - Comparison with a real-aperture radar 


The new Vespucci Bridge in Florence, Italy was inaugurated in 1957. It was 
not the reconstruction of the bridge destroyed by German mines during the 
WWII, but an original and innovative structure. Sixty years after its con- 
struction, the municipality commissioned an extensive monitoring campaign 
and possible restoration works. 

An interferometric real aperture (RA) radar has been used for testing 
and calibrate some sensors for a long-term monitoring. This case study has 
been exploited even for testing the interferometric MIMO radar. 

During the monitoring campaign of the Vespucci bridge the Ky-band 
RA (Ku-Ra) and the X-band MIMO (X-MIMO) radar have been positioned 
on the river bank close to a pillar (see map in Figure 4.36 and picture in 
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Figure 4.36). A Corner reflector (CR) has been installed on the lower deck 
(see pictures in Figure 4.36 and Figure 4.37(b)). 


(a) Picture of the two radars. (b) Picture of the view from the radars. 


Figure 4.37: Radar setup. 


Figure 4.38(a) shows the amplitude plot of the radar signal acquired by 
Ku-RA radar. The peak at about 25 m is the Corner reflector. The signal 
of the pillar in front of the radar is at about 47 m. 
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Figure 4.38(b) shows the image acquired by X-MIMO radar. The 2D 
imaging capability is demonstrated. The position of CR in the plan is cor- 
rectly detected. 
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(a) Amplitude plot acquired by Ku-RA (b) Radar image acquired by X-MIMO 
radar [128]. radar [128]. 
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Figure 4.38: Radar images. 


The main scope of this experimental test was to compare the displace- 
ment of CR detected by Ku-RA radar and X-MIMO radar. The sampling 
time of Ku-RA radar was 4.1 ms, while the sampling time of X-MIMO radar 
was about 30 s. 

Figure 4.39(a) shows the two plots in the same graph during the load 
operation of the bridge. Figure 4.39(b) shows the two plots during an- 
other measurement session, when the Ku-RA radar was operative only in 
the last 600 s. The agreement is rather good. 


4.4 Conclusions 


An interferometric compressive sensing MIMO has been designed and suc- 
cessfully tested as geotechnical equipment for testing bridges. 

The radar has been tested in controlled environment and the experi- 
mental results have been compared with the simulations. This test was 
performed by using a corner reflector on a micro-metric positioner stage. 

In-field test has been performed during a static test of a pedestrian bridge 
in Poggibonsi, Siena, Italy. The static load was located on the left side of the 
bridge and the radar was able to resolve different displacements of targets 
at the same range. 

The measurements exploited with the CS MIMO were also compared 
with a conventional real-aperture radar during the static test of Vespucci 
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(a) Comparison between Ku-RA amd X-MIMO 
radar [128]. Load operation. 
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(b) Comparison between Ku-RA amd X-MIMO 
radar [128]. 


Figure 4.39: Displacement of CR detected by RA radar (blue) and by MIMO 
radar (red) [128] 


bridge in Florence, Italy. The agreement between the results of two systems 
was good. 

The sampling time of the radar prototype is about 30 s, but this relatively 
long time is due to the specific VNA that operated as transceiver. It was an 
old model (HP8720D) not designed for fast operation. In effect, as this radar 
does not have mechanical moving parts, it could acquire much faster than 
the current GB synthetic aperture radar (SAR) based on the movement of 
a radar head along a mechanical linear guide. 

A fast interferometric radar can have single-tone integration time of 10 us 
with the number of frequency N; = 500. It means an acquisition time of 20 
ms using four antennas. A radar with 2-D imaging capability (provided by 
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the MIMO architecture) operating at this acquisition speed (50 Hz) opens 
exciting perspectives in the field of health monitoring and testing of large 
structure. 

For the sake of simplicity, the radar is provided with eight antennas. This 
configuration can give high side clutter with targets at more than 10°-20° 
with respect to direction of view. Nevertheless, a MIMO based on the same 
working principle could be used for monitoring targets with larger angular 
extension (like building or slopes) by increasing the number of antenna. 
Probably a CS MIMO with 8 x 8 antennas could cover the majority of 
applications. 
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Advanched in 3 dimensional GBSAR 


The current equipment exploits the movement of a radar head along a linear 
mechanical guide for synthesizing a large aperture [142]. Since the move- 
ment is along a single axis the obtained radar image does not have angular 
resolution in the plane orthogonal to the scan axis. In other words, if the 
radar head scans along the x-axis the radar image cannot have resolution 
in elevation angle. Figure 5.1 shows an example of resolution cell of the 
traditional systems. 

This is not a serious problem when the scenario is a slope, where the 
altitude (z-axis) can be reasonably considered an unambiguous function of 
the (x,y) position. Unfortunately there are cases where the geometry of the 
structure under test is much more complex [39], [79]. 

A trivial solution is to scan the radar head along the x-z plane (or on a 
spherical shell [143]). In this way it is possible to obtain a 3D image of the 
targets in the field of view. 

The problem is the measurement time can be very long. As an example, 
a fast GBSAR able to scan a line of 100 points in 5 minutes [87], will scan a 
plane of 100 lines in more than 8 hours. It means that the radar is able to 
detect only very slow movements. 


In this Chapter two radars able to provide a 3 dimensional image are 
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Figure 5.1: Resolution cell of traditional GBSAR systems. 


proposed. These radars are based on the technology proposed in the previous 
Chapters: the first uses the bistatic technique [144], [145]; the second model 
uses the CS method [146]. Recently, some other solutions for 3D imaging 
are presented for example [88] and [147]. 


5.1 GBSAR with 3D bistatic imaging capability 


A practical solution to achieve a 3D image could be a radar system able 
to acquire fast images without resolution in elevation and slow images with 
both azimuth and elevation resolution. The GBSAR presented in this section 
is a special bistatic GBSAR configuration that has this capability. 

The radar head moving along x-axis has a second transmitting channel 
linked to a third antenna through a RF cable. This third antenna is moving 
along a mechanical vertical axis. For each horizontal scan the radar system 
acquires a monostatic image of its field of view without elevation resolution. 
When the radar system completes the whole acquisition (along both x-axis 
and z-axis), it obtains a 3D image of its field of view. 

The sketch of the bistatic GBSAR system proposed is shown in Figure 
5.2. A linear monostatic GBSAR is moving along a horizontal rail (along x- 
axis) by acquiring interferometric 2-D images. The radar head has a second 
transmitting channel for operating in bistatic modality. By using this chan- 
nel, the radar transmits to a supplementary antenna that moves vertically 
along a mechanical axis. 
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The vertical axis is positioned at distance yo behind the horizontal axis 
in correspondence of its central point. The lower position of the antenna 
vertically moving is at zo over the x-axis. The radar head and the supple- 
mentary antenna are linked through a RF cable of length Lo. In this bistatic 
configuration the radar system is able to acquire data with spatial diversity 
in both direction (x and z) so it is able to obtain a 3D image of the targets 
in its field of view. 


ll 
N 
© --------------------» 


“o. 
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Figure 5.2: Working principle of the monostatic/bistatic GBSAR. 


The images can be focused in the 3D space in front of the radar, by taking 
into account of the path radar head - supplementary antenna - image point. 
The algorithm used is a modified version of the back-propagation algorithm 
described in (2.18). 

The advantage of this configuration with respect to other solutions [143], 
[148] is that for each horizontal scan the radar equipment is able to provide 
a monostatic image that can be used for detecting fast displacements. The 
complete 3D image is obtained only at the end of the vertical scan. 

Monostatic and bistatic images have different angular resolutions. As it 
is known, the azimuth resolution (Aym) of a monostatic GBSAR image is 
given by (2.14) 

À 
2Ly 


Apm = (5.1) 
with L, the horizontal span. 

With the aim to estimate the bistatic angular resolution, we have to take 
into account that transmitting from a point A and receiving from a point B is 
equivalent to transmit and receive from the median point C. With reference 
to Figure 5.3, the locus of median points is the dotted rectangle. 

Therefore, the bistatic images have angular resolution in azimuth (Agp) 
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Figure 5.3: Bistatic geometry. 


and elevation (A) given by 
À 
AnS oy. AS 


with W, and W, the two sides of the dotted rectangle in Figure 5.3. If the 
lower position of the supplementary antenna is close to the horizontal scan, 
W, = L,/2. Therefore 


(5.2) 


À À 
Ap x> — Ab, ~ — 5.3 
9° 7 = 7 (5.3) 
In other words, the azimuth resolution of the bistatic images is a factor 
2 worse than monostatic images. 


Radar prototype 


The GBSAR prototype we have assembled for testing the proposed working 
principle is shown in Figure 5.4. 

The vector network analyzer (VNA) HP8720D operated as transceiver 
providing a continuous wave stepped frequency signal (CWSF) in X-band 
with central frequency fe = 10 GHz and bandwidth B. Two RF cables linked 
the VNA to the front-end moving along the horizontal mechanical axis. Two 
single-pole double-throw (SPDT) switches provided a direct path between 
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Figure 5.4: Sketch the GBSAR prototype. 


the transmitter and the receiver in order to perform calibrated measure- 
ments. A third SPDT switched between the two TX antennas. 

For each single step along the horizontal mechanical rail the radar carries 
out the two acquisitions (monostatic and bistatic). A monostatic acquisition 
is complete after N, steps when the radar head has swept the horizontal axis, 
while the bistatic acquisition is complete only after N, horizontal scans, i.e. 
after N, x N; steps. The two mechanical axes spanned up to 1.50 m along 
x and up to 1.50 m along z. 


Experimental test 
Controlled environment 


In order to test the radar equipment in a controlled environment, the sensor 
has been installed on a garden. The map of the test field is shown in Figure 
5.5. 

A corner reflector (CR) of 0.4 m side was positioned on the tip of 3.0 m 
high wood pole. The pole was in front of the radar at 17.1 m distance. A 
second CR of 0.27 m side was positioned on a tripod (height from ground: 
1.2 m) on the right side of the radar. This CR could be moved along the 
view direction using a mechanical device with 0.1 mm nominal accuracy. 

The radar parameters are: f.= 10 GHz, B= 160 MHz, Ly= 1.5 m, 
L,= 1.02 m, N,= 150, N.= 102, yo= 0.63 m, zo= 0.19 m, Lo= 4.00 m. 

Figure 5.6 shows one of many monostatic images that have been obtained 
during the complete 3D bistatic acquisition. The acquisition time of each 
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Figure 5.5: Map of the test field. 


single monostatic image was 25 minutes, the complete 3D bistatic acquisition 
lasted about 42 hours. In the monostatic image are evident both the CR on 
the pole and the CR on the tripod, but obviously it is not possible to obtain 
their heights. 


Monostatic Image A.U.(dB) 


CR on tripod (1.2 m yy"? 


Figure 5.6: Monostatic radar image. 


During the bistatic acquisition the CR on the tripod was moved along 
the view direction using a micrometric positioner with 0.1 mm accuracy. 
Figure 5.7 shows the measured and the nominal displacement of the CR 
during the 102 acquisitions. The mean error was 0.13 mm. 

Finally, when the whole 3D bistatic acquisition was completed, it has 
been possible to focus in any arbitrary plane of the space in front of the 
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Figure 5.7: Measured (blue line with spots) and nominal displacement (red 
full line) of the CR on the tripod. 


radar. 


Figure 5.8 shows the bistatic image focused in the horizontal plane z = 
1.0 m, while Figure 5.9 shows the bistatic image focused in the horizontal 
plane z = 3.0 m. It is interesting to note that in the lower plane (z = 1.0 m) 
it is possible to see the CR on the tripod, the pole and the metallic fence, 
while in the upper plane (z = 3.0 m) plane only the CR on the tip of the 
pole is well visible, the CR on the tripod is not visible and the metallic fence 
gives only a light residual. 


CR on tripod (1.2 m ) -160 


0 
x (m) 


Figure 5.8: Bistatic radar image focused in the horizontal plane z = 1.0 m. 


Figure 5.10 shows the azimuth Point Spread Function of the corner re- 
flector on the tripod focused in the same plane both for monostatic and 
bistatic. It is possible to note that the resolution of bistatic is twice worst 
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Figure 5.9: Bistatic radar image focused in the horizontal plane z = 3.0 m. 


than monostatic. 
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Figure 5.10: Azimuth Point Spread Function of monostatic (red) and bistatic 
(blue) image. 


Figure 5.11 shows the bistatic image focused in the vertical plane at x= - 
1 m. It is evident the radar fence and the CR on the tip of the pole at 3.00 
m height. 


Realistic scenario 


With the aim of testing the equipment in a more realistic scenario, we pointed 
the radar forward a 7-storey building at about 140 m distance. 

The radar parameters we have set are central frequency 10 GHz, B= 80 MHz, 
L= 1.5 m, L,= 1.5 m, N,= 150, N,= 150, yo= 0.62 m, zo= 0.21 m, 
Lo= 4.00 m. 
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Figure 5.11: Bistatic radar image focused in the vertical plane x= -1.0 m. 


The obtained monostatic image is shown in Figure 5.12. The façade of the 
building is a clear oblique line. Obviously the monostatic (bi-dimensional) 
image cannot provide information about the height. 
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Figure 5.12: Monostatic image of the 7-storey building. 


Figure 5.13 shows the obtained radar image focuse in a horizontal plane 
z= 10.0 m. Also in this image the fagede are well visible as oblique line. 
Anyway, we can note that the resolution and the signal over noise are worst 
than Figure 5.12 

Bistatic radar image focused on the facede plane is shown in Figure 5.14. 
In this image the edges and the roof of the building are detectable. 

Finally, Figure 5.15 shows the radar image focused in the whole three- 
dimensional space in front of the radar. In order to remove the noise a 
threshold in amplitude equal to 1.7% of the maximum value was applied. 
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Figure 5.13: Bistatic radar image of the 7-storey building focused in the 
horizontal plane z= 10.0 m. 
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Figure 5.14: Bistatic radar image focused on the fagede plane of the 7-storey 
building. 


The fagade of the building is an evident target in the 3D space. 


5.2 GBSAR with 3D imaging capability using compressive sensing 
MIMO in elevation 


An other solution, for obtaining resolution in elevation, is to combine a 
multiple output multiple input (MIMO) radar [149] head able to provide 
vertical resolution with a horizontal scan. Indeed, this idea is inspired to 
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Figure 5.15: Bistatic 3D image. 


Synthetic Aperture Radar (SAR) tomography spaceborne [150] or airborne 
[151] that combine different passages for obtaining information in the third 
dimension. 

This hybrid approach has been tested in the same realistic measurement 
scenario of the preview system: a 7-storey building at about 140 m distance. 


Radar prototype 


The sketch of the radar system is shown in Figure 5.16. A step-frequency 
continuous wave (SFCW) transceiver transmits N} frequencies from fi to fa 
at step Af. The radar head is moving along a horizontal rail (the x-axis). 

The radar front-end was the same used in Chapter 4 for the interfero- 
metric CS MIMO. The pattern was the same of Chapter 4, though for this 
application the antennas where installed in vertical direction as shown in 
Figure 5.17. 

In this case the angular resolution (Ay) in the horizontal plane (azimuth) 
is obtained by mechanical scanning, while angular resolution (Av) in vertical 
axis (elevation) is obtained through the vertical disposition of the CS-MIMO. 
In both cases the angular resolution depends on the synthetic aperture: 


(5.4) 


(5.5) 
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Figure 5.16: Sketch of the radar system. 
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Figure 5.17: Configuration of antennas in MIMO radar head . 


For the radar system tested, L,=1.80 m, A=30 mm, N=40, so Ay = 
0.45° and AV’ = 2.86°. 
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Experimental test in a real measurement scenario 


The vector network analyzer (VNA) HP8720D operated as transceiver pro- 
viding a continuous wave stepped frequency signal (CWSF) in X-band. 
The radar parameters were: central frequency f.= 10 GHz, bandwidth 
B= 310 MHz, number of frequencies N;= 801, length of the horizontal 
scan L= 1.8 m, number of acquisition points along the x-scan N,= 180, 
transmitted power P= 26 dBm, gain of antennas G= 15 dB, time of inte- 
gration of each single tone r= 0.33 ms. The integration time was tint = 
Tx 16x Nf x Na = 769s. 

In order to test the radar equipment in a controlled environment, the 
radar has been installed on a garden of the University (see Figure 5.18). 
The aerial picture of the test field is shown in Figure 5.19. 


Figure 5.19: Aerial picture of the test site [146]. 


After a complete acquisition the measurement matrix Ek m can be fo- 
cused in the space in front of the radar. 

Figure 5.20(a) shows the obtained radar image focused on the horizon- 
tal plane at 10 m height and the radar image focused on a vertical plane 
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approximately on the façade of the farthest building (at about 140 m). 
The upper image in Figure 5.20(b) shows the radar image focused on a 


vertical plane approximately on the facade of the other building (at about 
100 m). 
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(a) Below: radar image focused on the (b) Below: radar image focused on the 
horizontal plane at 10 m height. Above: horizontal plane at 10 m height. Above: 
radar image focused on a vertical plane ap- radar image focused on a vertical plane 
proximately on the façade of the farthest approximately on the façade of the clos- 
building (at about 140 m) [146]. est building (at about 100 m) [146]. 


Figure 5.20: Radar images focused on the horizontal plane at 10 m height 
and radar images focused on a vertical plane approximately on the façade 
of the building. 


Finally, Figure 5.21 shows the 3D contour plot of the iso-surface using a 
suitable power threshold. The shape of the two buildings are clearly recog- 
nizable. 


5.3 Conclusions 


In this chapter two GBSAR prototypes with 3D imaging capability were 
proposed. 

The first GBSAR used the bistatic principle in order to obtain the eleva- 
tion resolution. This monostatic/bistatic GBSAR system has been demon- 
strated able to produce 2D images in short times (about a half hour) and 3D 
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Figure 5.21: The 3D contour plot of the power iso-surface with a suitable 
threshold [146]. 


images in long times (several hours). It is worth to note that this radar con- 
figuration has been obtained as an upgrade of a linear monostatic GBSAR. 
The bistatic 3D system shares the radar head and the horizontal mechan- 
ical axis with the 2D linear monostatic GBSAR. To improve the images 
resolution the central frequency can be increased(the prototype operated in 
X band, but the majority of GBSAR operates in Ky band). It is possible 
also to increasing the height of the vertical axis twice. In practice, a GB- 
SAR operating at 17 GHz with a 3 m vertical axis has vertical resolution 
four times better to the equipment presented in this thesis, and this can 
make the difference in many applications. 

The second GBSAR system used an hybrid working principle: synthetic 
aperture for azimuth resolution, and compressive sensing MIMO for elevation 
resolution. Also this system is able to provide a 3D image and the azimuth 
resolution can be improved by increasing the bandwidth and the scan length. 
On the contrary of azimuth resolution, elevation resolution does not depend 
on central frequency, but on the number of antennas used as shown in (5.5). 
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Conclusions 


The research activity, presented in this thesis, concerns design and develop- 
ment of some advanced Ground Based Radar Interferometric systems. These 
systems are related to the following issue: detection of displacement vector, 
Multiple Input Multiple Output (MIMO) and radars with 3D capability. 

The basic principle of this radar systems has been presented. All systems 
proposed have been successfully tested in controlled and real environment 
and, wherever possible, their performance have been compared with conven- 
tional systems. 


Radar systems for retrieving the displacement vector 


A GBSAR operating in monostatic and bistatic modality was presented. The 
sensor detects the first component of displacement as the conventional GBRI 
(monostatic) and an additional component through a transponder (bistatic). 

The radar has been tested with a basic transponder (two antennas and 
an amplifier) in a controlled environment and it was able to detect two 
components of displacement. The gain of the basic transponder was about 
50 dB. In this case the measurements were affected by some artefacts due to 
the signal liked from transponder during the monostitic measurements. 
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In order to improve the SNR of the bistatic measurements two enhanced 
transponders, with an hight gain, were presented. 

The first transponder made use of a frequency shifter in order to insulate 
the two antennas and achieved a gain of 91 dB. The artefacts did not affect 
the measurement with this transponder, but its implementation is rather 
difficult in practical cases. 

The other is a long cable transponder with a gain of 82 dB. This transpon- 
der appears to be more practical for most of the uses. 

Other application of bistatic technique was presented. A notable appli- 
cation is the monitoring of slander structures. In this application it was 
possible to measure the natural axis of a telecommunication tower for the 
first time with a GBRI. 

The issue of retrieving displacement vector is steel open. Indeed the 
GBRI presented in this thesis is able to detect only two components, the 
next step will be to generalize this theory in order to retrieve the whole 
displacement vector. 


Compressive Sensing MIMO Radar 


An interferometric compressive sensing MIMO has been designed and suc- 
cessfully tested as geotechnical equipment for bridges. The radar has been 
tested in controlled environment and the experimental results have been 
compared with the simulations. 

In-field test has been performed during a static test of bridge. The radar 
was able to resolve displacements of targets at the same range. 

The measurements exploited with the CS MIMO were also compared 
with a conventional real-aperture radar and the agreement was good. 

The sampling time of the radar prototype is about 30 s, but this relatively 
long time is due to the specific VNA that operated as transceiver. It was an 
old model (HP8720D) not designed for fast operation. In effect, as this radar 
does not have mechanical moving parts, it could acquire much faster than 
the current GB synthetic aperture radar (SAR) based on the movement of 
a radar head along a mechanical linear guide. 

A fast interferometric radar can have single-tone integration time of 10 
us with the number of frequency N; = 500. It means an acquisition time of 
20 ms using four antennas. A radar with 2-D imaging capability (provided 
by the MIMO architecture) operating at this acquisition speed (50 Hz) opens 
exciting perspectives in the field of health monitoring and testing of large 
structure. 

For the sake of simplicity, the radar is provided with eight antennas. This 
configuration can give high side clutter with targets at more than 10° — 20° 
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with respect to direction of view. Nevertheless, a MIMO based on the same 
working principle could be used for monitoring targets with larger angular 
extension (like building or slopes) by increasing the number of antenna. 
Probably a CS MIMO with 8 x 8 antennas could cover the majority of 
applications. 


Advanched in 3 dimensional GBSAR 


Two GBSAR prototypes with 3D imaging capability were proposed. 

The first GBSAR used the bistatic principle in order to obtain the eleva- 
tion resolution. This monostatic/bistatic GBSAR system was demonstrated 
able to produce 2D images in short times (about a half hour) and 3D images 
in long times (several hours). It is worth to note that this radar configura- 
tion has been obtained as an upgrade of a linear monostatic GBSAR. The 
bistatic 3D system shares the radar head and the horizontal mechanical axis 
with the 2D linear monostatic GBSAR. To improve the images resolution 
the central frequency can be increased(the prototype operated in X band, 
but the majority of GBSAR operates in Ky band). It is possible also to 
increasing the height of the vertical axis twice. In practice, a GBSAR oper- 
ating at 17 GHz with a vertical axis of 3 m has vertical resolution about four 
times better to the equipment presented in this thesis, and this can make 
the difference in many applications. 

The second GBSAR system used a hybrid working principle: synthetic 
aperture for azimuth resolution, and compressive sensing MIMO for elevation 
resolution. Also this system is able to provide a 3D image and the azimuth 
resolution can be improved by increasing the bandwidth and the scan length. 
On the contrary of azimuth resolution, elevation resolution does not depend 
on central frequency, but on the number of antennas used. 

The radar with 3D capability fascinates both the expert and the users. 
The solution presented in this thesis can be used for monitoring particular 
scenario while they are, in certain way, a retrofit of conventional GBSAR. 
This system, operating in higher (and wider) frequency could be used as 
high resolution 3D scanner. 


Industrialization 


The radar systems presented in this thesis are subject of three patents (prop- 
erty of the University of Florence). 

The patent related to bistatic radar was licensed by the University to a 
company, leader of GBRI sector, and probably it will be produce in short 
time. Regarding the patents on the CS-MIMO radar and ON the 3 di- 

125 


ADVANCED GROUND-BASED REAL AND SYNTHETIC APERTURE RADAR 


mensional GBSAR using monostatic/bistatic technique the University of 
Florence is currently negotiating for their licensing. 
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radars with 3D capability. 
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